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14.- Resumo/Notas 

Este projeto foi realtzado com os objetivos de desenvolver e 
tester o sistema de fiscalizagdo das operagoes de credito agricola, 
utilizando tecnicas de sensoriamento remote. A area piloto abrangeu os 
municipios de Aragatuba e Guararapes, localizados no Estado de Sdo Paulo. 
Utilizaram-se fotografias aereas, imagens LANDSAT e cartas topogrdficas. As 
fotografias aereas foram muito liteis para a delimitagdo das propriedades 
com financiamento do Banco do Brasil S.A. A percentagem de cana 

interpretada corretamente em imagens do LANDSAT, nas 85 propriedades 
analisadas, foi de 63,127a. A ocorrencia de condigoes atipicas como excesso 
de precipitagdo, cana bisada, florescida e tombada, verificadas no periodo 
de obtengdo de imagens, foram alguns dos fatores que contribuiram para o 
baixo desempenho da interpretagdo . Desenvolveu-se , entdo, uma abordagem 
alternative utilizando varies passagens do LANDSAT com suporte dos dados 
auxiliares que amentou a percentagem de acerto para 91,77%. As abordagens 
utilizadas (modelo agrometeoroldgico , indice vegetative e monitoramento da 
area de cana cortada) para estimative de rendimento, apresentaram 
resultados bastante encorajadores quando comparados com dados obtidos no 
campo. E dificil a deteegdo precise de areas menores que 25ha, exceto em 
condigoes favordveis de contexto. Os erros causados pela variagdo de escala 
de imagens podem ser considerados desvreziveis para pequenas areas 
qeogrdficas. 


15. Observagoes pesquisa financiada pelo banco do brasil s/a atravEs 

SEU FUNDO DE INCENTIVO A PESQUISA TECNICO-CIENTXFICA (FIPEC) . 
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ABSTRACT 


This project was carried out with the objectives of 
developing and testing remote sensing techniques for supporting the 
rural credit supervision system. The test area comprised the counties 
of Araoatuba and Guararapes, located in the State of Sdo Paulo. 
Aerial photographs, LANDSAT images and topographic charts were used. 
Aerial photograph:, were extremely useful for the cut lining of 
properties boundaries with financing of sugarcane plantations by the 
Banco do Brasil S.A. The percentage of correctly interpreted 
sugarcane on LANDSAT images, considering the 85 analised properties, 
was of 63.22%. The occurence of atypical conditions such as excessive 
raining, sugarcane in bloom, and wind damaged sugarcane and sugarcane 
not harvested due to planning failures verified during the period the 
images were obtained, were some of the contributing factors 
associated with a low interpretation performance. An alternative 
approach was developed using several LANDSAT overpasses and auxiliary 
field data, which resulted in 91.7? percent correct. The yield 
estimation approaches (agrometeorological model, vegetation index and 
monitoring of harvested sugarcane area) showed encouraging results 
when compared with field data. The accurate detection of areas 
smaller than 25 hi is difficult, except when having favorable 
conditions. The errors due to scale variations of the images can be 
neglected for small geographic areas. 


SUMARIO 


LISTA OE FIGURAS vii 

LISTA DE TABELAS ix 

CAPTTULO 1 - INTRQDUCAO 1 

CAPITULO 2 - MATERIAIS E METODOS 5 

2.1 - Materials 5 

2.1.1 - Descrigao geral da area piloto 5 

2.1.2 - Produtos de sensoriamento remote e material cartografico 6 

2. 1.2.1 - Fotografias aereas 6 

2. 1.2. 2 - Imagens do satelite LANDSAT 6 

2. 1.2. 3 - Documentagao cartografica 9 

2.1.3 - Equipamentos principals utilizados 10 

2. 1.3.1 - Zoom Transfer Scope (ZTS) 10 

2. 1.3. 2 - Ampliador fotografico Durst 10 

2. 1.3. 3 - Analisador de imagens multiespectrais (I-lOO) 10 

2.2 “ Metodologia para estimativa de area 11 

2.2.1 - Obtengao dos limites das propriedades agricolas 11 

2. 2. 1.1 - Confecgao de cartas tematicas 12 

2.2.2 - Obtengao e analise das fotografias aereas 13 

2. 2. 2.1 - Recobrimento aero fotografico 13 

2.2. 2.2 - Estabelecimento da legenda 14 

2. 2. 2. 3 - Estabelecimento de padroes 15 

2. 2. 2. 4 - Fotointerpretagao 16 

2. 2. 2. 5 - Revisao da fotointerpretagao e verificagao de campo ... 16 

2.2.3 - Obtengao e analise das imagens do LANDSAT 17 

2.2.3. 1 - Aquisigao das imagens do LANDSAT 17 

2. 2. 3. 2 - Estabelecimento de uma chave de interpretagao 17 

2. 2.3.3 - Interpretagao visual da cana-de-agucar (LANDSAT) 18 

2. 2. 3. 4 - Avaliagao de area 19 

2. 2. 3. 5 - Metodologia para verificar a exatidao do mapeamento da 

cana-de-agucar 19 


* - 


liECL’DING PAGE BLANK NOT FILMED 


Pag. 

2,2.4 - AnaHse dos dados do LANDSAT em fitas CCTs 21 

2.2.4. 1 - Procedimento para classificagao automatfca 21 

2. 2. 4. 2 - Estrategia para classificagao automatfca da area pil£ 

to 24 

2.3 - Metodologia para estimativa de produtfvtdade 2S 

2.3.1 - Modelo agrometeorologico 25 

2.3.2 - Obtengao de indices vegetativos para estimar rendimento 

da cana-de-agucar 27 

2.3.3 - Produtividade baseada la area de cana cortada 30 

CAPTTULO 3 - RESULTADOS E DISCUSSAQ 33 

3.1 - Dados obtidos atraves de fotografias aereas 33 

3.1.1 - Distribuigao espacial das propriedades agrTcolas 33 

3.1.2 - Precisao cartografica obtida no tragado dos limites das 

propriedades agrTcolas 39 

3.2 - Dados obtidos atraves de interpretagao de imagens do LANDSAT 42 

3.2.1 - Estimativa da area de cana-de-agucar cultivada por pro 

priedade agricola 42 

3.2.2 - Abordagem alternativa para identificagao de cana ao nTvel 

de propriedade agricola, atraves de imagens do LANDSAT .. 61 

3.2.3 - Tamanho minimo de area da cultura da cana-de-agucar d£ 

tetavel em imagens do LANDSAT 69 

3.2.4 - Avaliagao do uso da terra da area piloto 69 

3.3 - Classificagao automatica 71 

3.3.1 - Resultado da classificagao automatica 71 

3.4 - Modelos de produtividade 74 

3.4.1 - Produtividade atraves do modelo agrometeorologico 74 

3.4.2 - Produtividade atraves de indices vegetativos 77 

3.4.3 - Produtividade baseada na area de c^na cortada 78 

CAPTTULO 4 - CONCLUSDES 83 

4.1 - Ao nTvel municipal 83 

4.2 - Ao nTvel de propriedade 84 

CAPTTULO 5 - CONSIDERACOES FINAIS 85 


iv 


2hjL 

REFERENCIAS BIBLIOGRAFICAS 89 

APENDICE a - MAPAS DA DISTRIBUigAO DA CANA-DE-AgOCAR E ERROS ASSO 

CIADOS A SUA IDENTIFICAgAO 7 A.l 

APENDICE B - MAPA DE USD DA TERRA DA AREA PILOTO E.l 


LISTA DE FIGURAS 


Pag. 


1.1 - Area piloto que inclui os municTpios de Aragatuba e Guarar^ 

pes, situados na regiao noroeste do Estado de Sao Paulo ... 3 

2.1 - Esquana da cobertura aerofotografica da area de estudo 7 

2.2 - Fluxograma do procedimento 23 

3.1 - Regressao entre os valores de area de cana obtidos com d^ 

dos do LANDSAT e de aeronave 50 

3.2 - Caracterizagao de areas com problemas de cobertura de 

vens e perda de informagoes em imagers do LANDSAT 51 

3.3 - Regressao entre os valores da area de cana, obtidos com d^ 

dos do LANDSAT e dados de aeronave, referentes as 36 pro 
priedades agri colas que nao foram afetadas pel a falta de da 
dos do LANDSAT 7 53 

3.4 - Areas de cana cortada ate 26 de julho de 1983 nos munic^ 

pios de Aragatuba e Guararapes, obtidas da interpretagao de 
imagens do LANDSAT 79 

3.5 - Areas de cana cortada ate 27 de agosto de 1983 nos munic_T 

pios de Aragatuba e Guararapes, obtidas da interpretagao de 
imagens do LANDSAT 80 

3.6 - Areas de cana cortada ate 28 de setembro de 1983 nos municT 

pios de Aragatuba e Guararapes, obtidas da interpretagao de 
imagens do LANDSAT 81 


i 


iv i ■ 


- vtt - 






LISTA DE TABELAS 


Pag. 

2.1 - Dados do MSS do LANDSAT utilizados no projeto 8 

2.2 - Dados referentes a cobertura aerofotografica da area piloto 

compreendida pelos municTpios de Aragatuba e Guararapes ... 14 

2.3 - Calendario de plantio da cana-de-agucar nos municTpios de 

Aragatuba e Guararapes 20 

3.1 - Resultados obtidos referentes ao numero de propriedades 

agrTcolas, que cultivam a cana-de-agucar, financiadas pelo 
Banco do Brasil S.A., clentro da area piloto 33 

3.2 * Areas^das prcpriedades agrTcolas produtoras de cana,obtidas 

atraves de informagao verbal no campo e atraves de contagem 
de pontos sobre os mapas te-.aticos 1:250.000 34 

3.3 - Somatoria da area total, por destilarias, das propriedades 

agrTcolas produtoras de cana-de-agucar 38 

3.4 - Tabela demonstrative do erro devido a representagao grafica 

ao utilizer, para o desenho, pena de 0,1mm de diametro 40 

3.5 - Tabela demonstrativa do erro devido a representagao grafica 

ao utilizer, para o desenho, pena de 0,5mm de diametro 40 

3.6 - Tabela demonstrativa do erro devido a representagao grafica 

ao utilizer, para o desenho, pena de 0,5mm de diametro na 
escala de 1:50.000, reduzida para 1:250.000 41 

3.7 - Resultados comparatives de areas ocupadas com cana-de-agu 

car, por propriedade agrTcola, obtidos a partir de dados de 
aeronave e dados do LANDSAT 42 

3.8 - Resultados comparativos de areas ocupadas com cana-de-ag^ 

car, por propriedade agrTcola, obtidos por destilaria a 
partir de dados de aeronave e do LANDSAT 48 

3.9 - Resultado da interpretagao da cana-de-agucar a partir de 

imagens do LANDSAT e erros cometidos 54 

3.10 - Valores de cana interpretada corretamente e erros de incl^ 

sao e omissao por destilaria 59 

3.11 - Estimativa de area de cana, percentagem de cana interpret^ 

da corretamente e percentagem de erros de omissao e incl^ 
sao por destilaria 64 

3.12 - Proporgao da area ocupada^com cana-de-agucar, ao nTvel de 

propriedade, obtida atraves de imagens do LANDSAT e com 
trabalho simulado de campo, percentagem de cana interpret 
da corretamente e determinagao dos erros de omissao e iji 
clusao cometidos 65 

- ix - 


IMviiCEDlNG BLANK NOT FILMKD 



3.13 - Uso da terra da area piloto atraves de imagens do LANDSAT 71 

3.14 - Resultados do”PROCEDIMENTO-INPE 72 

3.15 - Valores medios de rendimento maximo (RM) e de rendimento 

real (RR), referentes a cana de ano e meio e a cana soca . 74 

3.16 “ Rendimentos medios em ton/ha de cana de ano e meio e de 

cana soca, em cada uma das cinco destil arias da area pilo 
to, obtidos atravis das informagoes das destilarias 75 

3.17 - Rendimentos de cana de ano e meio e de cana soca obtidos 

atraves do modelo agrometeorologico e das informagoes col^ 
tadas junto as destilarias da area piloto 76 

'3.18 - Coeficientes de correlagao (r) entre os valores de rendj_ 
mento da cana-de-agucar, dados espectrais e indices vege 
tativos 77 

3.19 - Rendimento medio de cana das destilarias de Aragatuba e 

Guararapes 82 


CAPTTULO 1 


introducao 


A agricultura tem sido considerada como urn dos setores 
prioritarios do Governo que participa significativamente da economia 
do Pais. Considerando a extensao tei^ritorial do Brasil, aliada a 
necessidade de aumento de producao devido ao crescimento 
populacional , verif1ca-se utna acentuada tendencia de expansao da 
fronteira agricola. Assim sendo, torna-se necessaria a adocao de 
tecnicas racionais e eficientes para o controle e monitoramento dos 
processos de producao agrTcola. 

0 Banco do Brasil S.A., atraves da operajao de credit© 
agrTcola, tem incentivado e possibil itado aos agricultores aumentar 
a area plantada e a produtividade, gerando, consequentemente, uma 
maior quantidade de produtos. Para garantir a efetividade do credito 
agrTcola, torna-se de fundamental importancia o estabelecimento de urn 
sistema eficiente de fiscal izacao dos recursos aplicados para o 
foment© da Agricultura. A fiscal izacao, baseada em metodos 
convencionais, i lirritada devido ao dinamismo da atividade agrTcola e 
a rapidez de expansao das areas cultivadas, distribuindo-se em 
grandes regioes do PaTs. Alem disso, metodos convencionais sao, em 
geral, bastante deficientes na estimativa de areas no campo- Urn 
sistema baseado em sensoriamento reuioto apresenta caracterTsticas que 
poderao ser de grande utilidade no process© de fiscal izacao. 

0 sensoriamento remote e uma tecnica que consiste na 
obtencao de dados da superfTcie da Terra, por meio de instrumentos 
colocados, principalmente, a bordo de aeronave e de satelite. 0 
lancamento de satelites de observacao da • Terra e o crescente 
desenvolvimento das tecnicas de anilise tim permitido coletar dados 
da superfTcie terrestre, possTveis de ser tratados de modo rapido, 
eficiente e economico, gerando informacoes altamente estrategicas e 
favorecendo o controle dos recursos terrestres. Pel as suas 
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caracteristicas de recobrir periodicamente grandes regioes,; o 
sensorimento remote ao nivel orbital permite que se estabeleca uma 

f 

divers idade de programas de monitoramento dos recursos da Terra de 
forma dinamica e rapida. 

Este projeto visou desenvolver e testar urn siscema de 
monitoramento de cultura ao nTvel de municTpio e de propriedades 
agrTcolas, utilizando tecnicas de sensoriamento remote, com enfase na 
utilizacao de dados coletados por satelites de observagao da Terra. 

Como area piloto, foram escolhidos os municTpios de 
Aracatuba e Guararapes, na regiao noroeste do Estado de Sao Paulo 
(Figura 1.1), ocupando uma area de 3.583km^, regiao de grande 
tradicao de exploracao pecuaria, mas que atualmente apresenta 
acentuada e rapida expansao agricola, principalmente da cultura da 
cana-de-acucar. 
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Fig. 1.1 - Area piloto que inclui os municTpios de Aragatuba e Guarara 
pes, situados na regiao noroeste do Estado de Sao Paulo. 
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CAPITULO 2 


MATERIAIS E METODCS 


2.1 - MATERIAIS 

2.1.1 - descricao geral da Area piloto 


A area piloto compreende os municTpios de Aracatuba e 
Guararapes. Esta loca.izada entre as latitudes 20®45 e 21°30 sul e 
longitudes 50 15 e 51 00 oeste (Figura 1.1). Situada a noroeste do 
Estado de Sao Paulo, ocupa uma superfTcie de 358.300ha, abrangendo 
areas agrTcolas das destilarias Cruzalcool, Alcoazul, Aralco e 
Destivale, local izadas no municTcpio de Aracatuba e da destilarla 
Unialco, local izada no municTpio de Guararapes. 

Em termos de caracterTsticas naturals, a area escolhida, 
insere-se na ProvTncia Fisiografica denominada Planalto Ocidental 
Paulista, que apresenta relevo com interfluvios amplos e suavizados, 
constituindo-se em divisores hidrograficos dos rios Tiete e AguapeT. 
Em termos de solos, predominam os solos do grupo Latossolo Vermelho 
Escuro Fase arenosa e os solos Podzolizados de Lins e MarTlia (Brasil, 
1960). Sobre as caracterTsticas de clima, predomina na area o tipo 
climatico Awa de Kdppen, subtropical, seco no inverno 
(junho-julho-agosto) e verao chuvoso (novembro-dezembrc-janeiro) , com 
temperatura media do mes mais quente superior a 22°C. A precipitacao 
anual esta er.tr? 1000 e 1500mm e a temperatura media anual entre 20 e 
22°C (Setzer, 1966). 

Do ponto de vist? agrTcola, a pecuaria e a exploracao 
dominante na regiao, a qual ocupa vastas areas de pastagens. A cultura 
da cana-de-acucar nao e expressiva na area de estudo, cuja 
implantagao, para fins de exploracao industrial, teve inTcio no ano de 
1980. 
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2.1.2 - PRQDUTOS DE SENSORIAMENTO REMOTO E MATERIAL CARTOGRAFICO 

% 

2. 1.2.1 - FOTOGRAFIAS AFREAS 


Foram utilizadas fotografias aereas pancromaticas 
verticals, resultantes das coberturas aerofotograficas realizadas, em 
1972, pela FAB 1/69 GAV, para o Institute Brasileiro do Cafe (IBC) e, 
em 1978, pela Terrafoto S.A. para a Companhia de Eletrificacao de Sao 
Paulo (CESP). Estas fotografias foram utilizadas, em campo, para 
delimitacao das propriedades agrTcolas levantadas (Figura 2.1). 

Para a interpretacao da cana-de-acucar foram analisadas 
transparincias de fotografias aereas infravermelhas coloridas 
verticais na escala 1:30.000, resultantes da cobertura 
aerofotografica realizada especialmente para esse projeto, pelo 
Institute de Pesquisas Espaciais (INPE), no periodo de 13 a 23 de 
junno de 1°83, cobrindo integralmente a area dos dois municipios 
(Figura 1.1). 

2. 1.2. 2 - IMAGENS DO SATELITE LANDSAT 

0 satelite da sirie LANDSAT possui urn imageador 
multiespectral (MSS-Multispectral Scanner Subsystem), o qual obtem 
simultaneamente imagens da superfTcie da Terra em quatro canais 
espectrais: canal 4: 0,5 a 0,6ym (verde-laranja) i canal 5; 0,6 a 0,7ym 
(laranja-vermelho) ; canal 6: 0,7 a 0,8ym (vermelho-infravermelho 
proximo); canal 7: 0,8ym a l,1ym (infravermelho proximo). 

As informacoes obtidas por esses sensores sao enviadas 
para a Terra, na forma de sinais eletromagneticos, que sao recebidos e 
gravados por estacoes de recepcao. Posteriormente, sao' convertidos em 
imagens fotograficas e fitas digitais compatTveis com computador. 
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Fig. 2.1 - Esquema de cobertura aerofotografica da area de estudo. 





No presente projeto foram estudadas as imagens obtidas 
pelo imageador multiespectral (MSS), processadas nos laboratories do 
Institute de Pesquisas Espaciais (INPE), identificadas pela orbita 222 
do satelite LANDSAT nos pontos 74 e 75, analisadas sob duas formas de 
produto, conforme mostra a Tabela 2.1. 

TABELA 2.1 




















2. 1.2. 3 - DOCUMENTACSO CARTOGRAFICA 


• Para a identificacao de acidentes geograficos nas 
fotografias aereas e imagens do LANDSAT, foram utilizadas como base as 
cartas topograficas nas escalas 1:50.000 e 1:250.000, conforme relacao 
abalxo: 


Cartas Topograficas, escala 1:50.000: 


Sud Mennucci 

- 

1967, 

I6G/SP 

Auriflama 

- 

1967, 

I6G/SP 

Bacuri 

- 

1967, 

IGG/SP 

Aracangua 

- 

1967, 

IGG/SP 

Major Prado 

- 

1967, 

IGG/SP 

Gastao Vidigal 

- 

1967, 

IGG/SP 

Valparaiso 

- 

1^67, 

IGG/SP 

Ribeiro do Vale 

- 

1967, 

IGG/SP 

Aracatuba 

- 

1967, 

IGG/SP 

Bento de Abreu 

- 

1967, 

IGG/SP 

Guararapes 

- 

1957, 

IGG/SP 

Birigui 

- 

1967, 

IGG/SP 

Salmourao 

- 

1974, 

FIBGE 

Rinopol is 

- 

1974, 

FIBGE 


- Cartas Topograficas, escala 1:250.000; 

Votuporanga - 1976, FIBGE 

Aracatuba - 1976, FIBGE 

Obteve-se, tamben, para servir como fundo da 
interpretacao visual da cana-de-acucar a partir de imagem de satelite, 
uma base planimetrica extraida de uma imagem de satelite, que foi 
geometricamente corrigida, denominada Precision. 


2.1.3 - EQUIPAMENTOS PRINCIPAIS IfTlLIZADOS 


2.1.3. 1 - ZOOM TRANSFER SCOPE (ZTS) 


Para a transferincia dos limites das propriedades, 
contidas nas fotografias aereas em 1:20.000 e 1:25.000 sobre as 
cartas topograficas em 1:50.000, foi utilizado o ZTS (Zoom Transfer 
Scope) que permite urn ajuste de escala atraves da superposicao de 
acidentes geograficos. 

2. 1.3. 2 - AMPLIADOR FOTOGRAFICO DURST 


Para promover urn ajuste grafico dos mapas tematicos em 
1:50.000 para a escala de trabalho de 1:250.000 foi utilizado o 
ampliador fotografico Durst do laboratorio fotografico do INPE. 

2. 1.3. 3 - ANALISADQR DE IMAGENS MULTIESPECTRAIS (1-100) 

0 1-100 e urn sistema interativo, conversacional , 
projetado para extrair automat icamente informacoes multiespectrais de 
imagens terrestres obtidas remotamente por satelites ou avioes. 
Baseia-se no princTpio de que todos os objetos possuem 
caracterTsticas espectrais unicas que sao tambem denominadas de 
assinaturas, isto e, cada objeto reflete ou emite energia 
eletromagnetica de uma maneira peculiar nas diversas bandas do 
espectro. 


0 1-100 possui urn "hardware" especial que opera sob o 
controle de urn minicomputador PDP-11/45 que controla tambem os 
algoritmos (software) implementados, os quais permitem identificar 
areas semelhantes na imagem sob estudo, comparando assinaturas 
previamente conhecidas com as caracterTsticas espectrais de cada 
ponto da imagem, associando cores falsas as areas da imagem que 
possuem as mesmas caracterTsticas espectrais (temas). Estes temas sao 
entao visual izados em urn monitor de televisao, individualmente ou 
simultaneamente com a imagem original. 


2.2 - METODOLOGIA PARA ESTIMATIVA DE AREA 


2.2.1 - OBTENCAO DOS LIMITES DAS PROPRIEDADES AGRiCOLAS 

Para a obtencao em campo dos limites das propriedades, 
realizou-se preliminarmente uin teste para definir quais os documentos 
cartograficos ou fotograficos que poderiam ser utilizados. Os 
documentos a disposiqao foran os seguintes: 

a) cartas topograficas, escala 1:250.000; 

b) cartas topograficas, escala 1:50.000; 

c) imagens MSS do LANDSAT, escala 1:250.000; 

d) imagens RBV do LANDSAT, escala 1:100.000; 

e) fotografias aereas pancromaticas, escala 1:20.000, voo de 

1978; 

f) fotografias aereas pancromaticas, escala 1:25.000, voo de 

1972. 


Definido o material que melhor se prestava ao objetivo 
de delimitacao das proprieoades (no caso o material dos itens "e" e 
"f"), procedeu-se, sistematicamente, ao levantamento das propriedades 
adotando os seguintes passos: 

1) Com auxTlio de urn fotoTndice, verif icava-se em qual linha de 
voo estava local izada a propriedade, cujo proprietario 
deveria ser entrevistado. A seguir, selecionavam-se as 
aerofotograf ias que cobriam a propriedade a ser delimitada; 

2) Com as aerofotografias selecionadas, compunha-se urn pequeno 
mosaico e mostrava-se ao entrevistado detalhes da regiao, 
registrados sobre as aerofotografias, a fim de que ele 
pudesse local izar a sue propriedade nas aerofotos; 
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3) A seguir, com a orientacao do proprietario, tracava-se sobre 
as fotografias aereas, com lapis do tipo cera, os limites da 
propriedade;. 

4) Tracados os limites da propriedade, marcava-se sobre a 
fotografia aerea, ou conjunto de fotografias aereas, o numero 
de codigo de identificacao daquela propriedade demarcada; e, 

5) Finalmente, sobre uma folha de papel , fazia-se a descricao da 
propriedade, com os seguintes dados: 

- numero de codigo da propriedade, 

- nome do proprietario, 

- area da propriedade, 

- numero das fotos aereas nas qua is a propriedade estava 
contida. 

2.2. 1.1 - CONFECCAO de CARTAS TEMATICAS 


0 lancamento dos limites das propriedades agrTcolas 
contidas nas fotografias aereas referidas no item 2.2.1 (e,f), foi 
feito sobre bases cartograficas na escala de 1:50.000, utilizando uma 
camara clara ZTS (Zoom Transfer Scope). 

A interpretacao das fotografias aereas referidas na 
Secao 2. 1.2.1, paragrafo 2, foi reduzida fotograficamente e ajustada 
sobre as bases cartograficas na escala de 1:50.000. 

0 desenho das cartas tematicas com os limites das 
propriedades agrTcolas e a cultura da cana-de-acucar, foi obtido por 
compilaqao direta das folhas topograficas na escala de 1:50.000, 
elaboradas pelo Institute Geologico e Geografico e Fundaqao Institute 
Brasileiro de Geografia e EstatTstica (IGG/SP e FIBGE). 


y 


- 13 - 


Confeccionadas as cartas tematicas na escala de 

1:50.000, procedeu-se a sua reducao para a escala de 1:250.000 a fim 
de permitir a comprovacao entre os dados interpretados pelo LANDSAT e 
os das fotografias aereas (verdade terrestre). 

Tais produtos podem ser vistos no atlas que segue anexo 
a este relat5rio. 

2.2.2 - OBTENCAO E ANALISE DAS FOTOGRAFIAS AEREAS 
2.2.2. 1 - RECOBRIMENTO AEROFOTOGRAFICO 

A cobertura aerofotografica foi executada com os 

seguintes objetivos: 

1) auxiliar no estabelecimento de padroes fotograficos com o 
objetivo de estabelecer uma chave de interpretacao baseada em 
areas de treinamento representativas dos temas a serem 
identificados; 

2) servir corao verdade terrestre pare medir a exatidao de 
identificacao e avaliacao de area de temas interpretados em 
image ns LANDSAT. 

Para a realizacao desse voo aerofotografico foi 
utilizada a aeronave Bandeirante do INPE, no perTodo de 13 a 23 de 
junho de 1°^3. A Tabela 2.2 mostra os dados tecnicos da cobertura 
aerofotogr if ica da area de estudo. 
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DADOS REFERENTES A COBERTJRA AEROPOTOGRAFICA DA AREA PHOTO 
COMPREENDIDA PELOS MUNICTPIOS DE ARACATUBA E GUARARAPES 


cAmara 

distancia 

FOCAL (mm) 

FILME 

(KODAK) 

FILTROS 

RECOBRIMENTO 

NO DE 
FAIXAS 

NQ DE rOTOS 

ALTURA 

(PES) 

LONG. 

LAT. 

TOTAL 

OTEIS 

WILD 

RC-10 

152 

Aero- 
chrome 
I R- 2443 

W500 + 
Av.1,4 

60^ 

30% 

21 

466 

389 

15000 


2. 2. 2. 2. - ESTABELECIMENTO DA LEGENDA 


Para os primeiros contatos com os padroes de imagem da 
area piloto, foram selecionadas quatro areas de treinamento de 
10km X 10km, dentro da area de estudo, distribuidas segundo criterio 
de diferentes concentracoes da cultura da cana, conforme pode ser 
visto na Figura 1.1. As fotografias aereas correspondentes a estas 
areas de treinamento foram interpretadas pelos elementos da equipe. 
Os resultados foram discutidos e a area piloto foi visitada para dar 
inTcio a correlaqao entre os padroes das aerofotos e o terreno. 

Como resultado desta fase, foi estabelecida a seguinte 

legenda: 

C - Cana adulta (Safra 1983), 

CC - Cana corlada (Safra 1983), 

CN - Cana Nova (Plantio 1983). 

Nos mapas do atlas anexo, que contem a distribuipao 
espacial da cultura canavieira, na area piloto, a legenda acima foi 
assim simpl ificada: 


I 

I 

I 


• Vi 
\ 

1 

m- 

T 
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C + CC - Cana safra 1983, 

CN - Cana plantio 1983. 

2.2. 2.3 - ESTABELECIMENTO DE PADROES 


Foram estabelecidos os padroes e as normas para 
orientacao durante os trabalhos de interpretagao final, com base na 
experiencia adquirida nas fases anteriores. 

a) CAM ADULTA - Em estigio avangado de desenvolvimento (a ser 
cortada na safra em andamento), na epjca da tomada da 
fotografia aerea. Dentro da classe da cana adulta incluem-se 
cana de ano e meio, cana de ano e cana soca, aqui definidas 
atraves de padroes de tonal idade e textura: 

- Tonalidade: Magenta mais ou menos intensa, dependendo da 
manor ou menor percentagem de cobertura do solo pel a 
planta. No caso da ocorrencia de condigoes anormais, 
devido a fatores ambientais e agronomicos, a tonalidade se 
apresentava magenta, mosqueada de manchas brancas. A 
variedade plantada tambem influ iu na variagao tonal. 

- Textura: Aveludada de fina a grosseira, dependendo da 
ocorrencia de cana de ano, cana soca e cana de ano e meio. 

b) CANA CORTADA - Areas col hi das durante a safra em andamento 
na epoca da cobertura aerofotograf ica : 

- Tonalidade: Verde com faixas brancas causadas pela 

presenga de alinhamentos em paralelo das leiras de palhada 
que acompanham o sentido do corte. 

- Textura: Fina e descontTnua (aspecto penteado). 

c) CAM NOVA - Cana de ano e meio em inicio de desenvolvimento, 
a qual nao sera cortada na safra em andamento. 
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- Tonalidade: Variivel , conforme 0 estagio da cana e o tipo 

do solo. Em geral, urn padrao de verde intense para solos 
« 

de coloracao vermelho dominante, e urn verde claro para 
solos de coloragao clara. 

- Textura: Fina e continua. 

2. 2. 2. 4 - FOTOINTERPRETACAO 


A fotointerpretacao foi desenvolvida atraves do uso 
de mesas de luz e lupas monoculares com aumento de 8 vezes, as quais 
propiciam a observacao dos aspectos de tonalidade e textura, 
criterios altamente relevantes para a identificagao da cultura da 
cana-de-agucar em transparencias coloridas de fotografias aereas 
infravermelhas. 


Os trabalhos da fotointerpretagao propriamente ditos 
foram antecedidos pela marcagao de pontos homologos nos pares de 
fotografias aereas contiguas, com o objetivo de minimizar as 
distorgoes que ocorrem em suas bordas e otimizar o trabalho de 
fotointerpretagao, evitando que a mesma area fosse interpretada mais 
de uma vez. 

2. 2. 2. 5 - REVISAO DA FOTOINTERPRETACAO E VERIFICACAO DE CAMPO 


Todos os pontos com padroes duvidosos, resultantes da 
fotointerpretagao realizada individualmente, eram submetidos i 
apreciagao, para decisao, dos demais membros da equipe. 

Concluidos os trabalhos de fotointerpretagao, foram 
ioentificados sobre algumas fotografias pontos com padroes 
indefinidos e pontos com padroes bem definidos a serem visitados no 
campo, a fim de dirimir duvidas e confirmar decisoes tomadas. 
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2.?. 3 - OBTENQAO E ANALISE DAS IMAGENS DO LANDSAT 
2.2.3. 1 - AQUISICAO DAS IMAGENS DO LANDSAT 


A aquisipao das imagens do LANDSAT foi feita levando em 
consideracao dois criterios: 

1) qualidade das imagens, em termos de cobertura de nuvens; e 

2) perTodo de aquisicao das imagens, em funcao do calendar io 
agricola da cana-de-acucar adotado pelos agricultores da 
regiao. 

2.2.3. 2 - ESTABELECIMENTO DE UMA CHAVE DE INTERPRETACAO 


0 plantio da cana-de-acucar no Estado de Sao Paulo, 
segundo Koffler et alii (1979), e concentrado en duas epocas. 0 
plantio feito em outubro, que permite a colheita de cana com 
aproximadamente 12 meses (a partir de setembro do ano seguinte), 
sendo conhecida como cana de ano, e 0 plantio nos meses de Janeiro e 
fevereiro, que permite a colheita de cana com 18 meses (a partir de 
maio do ano seguinte), conhecida como aana de ano e meio. 

Com base neste esquema de plantio e corte da 
cana-de-acucar, foi utilizada uma chave de interpretagao desenvolvida 
e testada por Mendonca et alii (1981) para interpretar 
a cana-de-acucar da area deste Projeto. 0 princTpio basico desta 
chave e a utilizagao combinada das caracteristicas espectrais e 
temporais das imagens do LANDSAT para identificar a cana em seus 
diferentes escagios de desenvolvimento. 

Deve-se ressaltar, entretanto, que esta chave de 
interpretacao foi desenvolvida para situaqoes em que a cana apresenta 
urn desenvolvimento normal, sem a ocorrincia de problemas que venham 
alterar a sua resposta espectral caracterTstica. 


Durante o ano de 1983 ocorreram na regiao de estudo 
problemas agronornicos e meteorologicos tais como: cana tombada, cana 
praguejada e com inflorescencia, cana bisada, alem de ocorrerem 
Indices de precipitacao acima do normal para aquele perTodo de 
analise. Esta situacao provocou mudancas no comportamento espectral 
da cana-de>acucar, o que impediu a utilizapao da chave de 
interpretagao acima citada. 

Para contornar este problema foi desenvolvida outra 
chave que levou em consideragao a caracterizagao das areas preparadas 
para o plantio da cana e datas de plantio especTficas para a area de 
estudo. Esta alternativa implicou um maior numero de imageis. 

2.2.3. 3 - INTERPRETACAO VISUAL DA CANA-DE-AQOCAR (LANDSAT) 


Os trabalhos de identificagao da cana-de-agucar , a 
partir de imagens do LANDSAT, sao referentes a 85 propriedades 
agricolas, correspondentes as propriedades fornecedoras de cana para 
as cinco destilarias da area de estudo (Tabela 3.7), e cujas areas de 
plantio de cana sao financiadas pelo Banco do Brasil. 

A identificagao e delimitagao das areas de solo 
preparado foram posslveis devido ao grande contraste existente entre 
0 solo exposto e os alvos adjacentes. Alim disto as areas de solo 
preparado apresentam formas geoni-tricas bem definidas. 

Para esta regiao, as areas de solo preparado se 
apresentavam, no canal 5, com nlveis de cinza claro. As datas de 
passagens das imagens utilizadas foram as seguintes: 


,\no de 1980; 
Ano de 1981 ; 
Ano de 1982; 
Ano de 1983: 


16/jan.; 

14/mar. ; 

14/jan., 18/mar., 23/abr., 29/maio e 22/jul.; 
16/fev., 07/maio e 10/jul. 
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Para decidir quais, dentro das areas de solo 
preparado, eratn aquelas provavelmente destinadas ao plantio da 
cana-de-acucar, foram utilizadas imagens do LANDSAT, no canal 7, 
referentes as seguintes passagens: 

Ano de 1982; 29/maio e 22/jul.; 

Ano de 1983; 07/naio e 10 /'j1. 

A escolha destas datas foi baseada em informacoes 
. obtidas em campo sobre o calendario de plantio da cultura naquela 
regiao, desde o ano de sua implantagao ate o ano de 1983, mostrado na 
Tabela 2.3. 

2. 2. 3. 4 - AVALIACAO DE AREA 


0 metodo utilizado para avaliacao de area de c:''sa foi o 
de coritagem de pontos de uma malha transparente sobreposta aos mapas 
tematicos. 


Para o caso de mapas tematicos, na escala de 1:250.000, 
obtidos de imagens do LANDSAT, foi usada grade de 1mm x 1mm; enquanto 
para os mapas de 1:50.000, gerados a partir de fotografias aereas, 
utilizou-se grade de 0,5cm x 0,5cm. Cada ponto corresponde a 6,25 ha 
em ambas as grades. 

2.2.3. 5 - METODOLQGIA PARA VERIFICAR A EXATIDAO DO MAPEAI^ENTO DA 

cana-de-acdcar 


Para verificar a exatidao do mapeamento da 
cana-de-acucar, utilizando imagens do LANDSAT, adotou-se . 
primeiramente os seguintes criterios preestabelecidos; 

1) que 0 risco do consumidor em aceitar urn mapa com uma exatidao 
inferior a 85% seja de 5%. 


CALENDRRIO OE PLANTIO DE CANA-DE-ACDCAR NOS MUNIClPIOS DE ARACATUBA E GUARARAPES 


i 


ORIGINAL PAGE 19 
OF POOR QUALITY 






risco do consumidor (0,05) = z 111 0.85*^“^ (1-0.85)^; (2.1) 

Y=0 Y! (N-Y)! 


2) que 0 risco do produtor em rejeitar um mapa com uma exatidao 
peri or a 95% seja tambem de 5%: 


risco do produtor (0,05) = z (1-0.95)^; (2.2) 

Y = X+ 1 Y! (N-Y)! 


onde: 


N = numero de pontos amostrados, 

Y = numero de pontos erroneamerte classificados , 

X = numero maximo de pontos erroneamente classificados para sati^ 
fazer as Equa?oes 2.1 e 2.2. 

De acordo com estes criterios, foi dimensionada uma amos 
tra minima de 93 pontos a serem checados. A aceitagao do mapa, dein 
tro dos criterios estabelecidos , estaria condicionada ao fato de nao 
mais do que 8 pontos, dentro dos 93 selecionados , serem incorretos. 

2.2.4 - ANALISE DOS DADOS DO LANDSAT EM FITAS CCTs 

2. 2. 4.1 - PROCEDIMENTO PARA CLASSIFICACAO AUTOMATICA 


Com 0 objetivo de desenvolver um sistema integrado de pr£ 
visao de safras, foi elaborado no INPE um procedimento interativo que 
envoi ve analista e computador na determinagao da proporgao de culturas 
agrTcolas em segmentos de imagers de satellte .(Caniara et alii, 1982). 
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A metodologia proposta para a estimativa de area 
atraves do uso deste procedimento i a de amostragem estratif icada e 
envolve os seguintes passsos: 

1) zoneamento da regiao de estudo de acordo com as 
caracteristicas edaficas, climaticas e agricolas, 

2) estratif icagao das zonas homogeneas de acordo com a 
densidade da cultura em estudo; 

3) selecao aleatoria de unidades amostrais, denominadas 
segmentos, em cada estrato homogeneo; 

4) determinacao da proporcao da cultura em cada segmento; 

5) estimativa da area da cultura por estrato atraves da 
proporcao estimada nos segmentos; 

6) as estimativas dos estratos sao agregadas e obtem-se a 
estimativa para as zonas homogeneas que, por sua vez, sao 
agregadas para obter a estimativa da cultura na regiao de 
estudo. 

0 procedimento utiMzado no presente trabalho envolve 
a participacao do analista (agrcnomo fotointerprete) e do computador, 
que trabalham de forma interativa, com quatro fases distintas: 

a) treinamento e rotulacao; 

b) classificacao; 

c) estimativa de proporcao; 

d) aval iacao e aprovacao, 


conforme mostra a Figura 2.2, detalhadas a seguir. 
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Fig. 2.2 - Fluxograma de procedimento. 


a) Treinamento e rotulagao 

Nesta fase o analista rotula manualmente pontos 
prialocados aleatoriamente sobre a imagem com dois objetivos; fornecer 
os parametros necessarios para a classificapao; dar subsTdios para 
testar a estimativa da proporgao e corrigir erros de tendencia do 
classif icador na estimativa dessa proporgao. 

b) Classificagao 

A classificagao separa cada ponto do segmento em duas 
classes: cultura e nao cultura, de acordo com a regra de maxima 

verossimilhanga. 

c) Estimativa de proporgao 

0 resultado da classif icagao e usado para a estimativa 
da proporgao da cultura no segmento. Nesta fase tambem e feita uma 
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corregao das eventuais tendencias ocorridas na classificacao que 
utiliza os subsidies fornecidos na fase de treinamento e rotulacao. 

u) Avaliacao e aprovacao 

0 resultado da estimativa da proporcao e submetido ao 
analista para aprovacao e, a seu criterio, os passes do procedimento 
poderao ser refeitos, ou o segmento estara pronto para fornecer a 
proporcao da cultura plantada. 

As principals caracterTsticas do procedimento sao: 

- minimizacao do uso de informacoes de campo; 

- capacidade de verificar os resultados da classificacao 
automatica com a rotulacao feita pelo analista. 

2.2. 4. 2 - ESTRATEGIA PARA CLASSIFICACAO AUTOMATICA DA AREA PHOTO 


Na area piloto foram alocados seis segmentos, 
denominados segmentos testes, nas proximidades das destil arias da 
regiao (Figura 1.1), 

Inicialmente foram selecionadas quatro passagens do 
satelite LANDSAT, cujos dados contidos em fitas CCTs (fitas 
compativeis com o computador) foram digitalmente sobrepostos e 
transfer idos para uma unica fita de computador. 

Para cada urn dos segmentos foram processadas as fases de 
treinamento e rotulacao, classificacao, estimativa da proporcao, 
avaliacao e aprovacao. 

As passagens utilizadas para estimativa da proporcao da 
cultura nos segmentos foram 16 de fevereiro de 1983, 07 de maio de 
1983, 10 de junho de 1983 e 27 de agosto de 1983, denominadas 
passagens 1, 2, 3 e 4, respectivamente. 
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Somente para os segmentos 2 e 4 estavam disponTveis as 
quatro passagens. Para o segmento 1 nao se dispunha da passagem 3 e os 
segmentos 3 e 5 estavam com cobertura de nuvens na 2? passagem. 0 
segmento 6 nao foi utilizado para c^assificagao pois nao apresentava 
nenhuma area de cana. 

0 procedimento visa obter a estimativa da area plantada 
com culturas agricolas em grandes areas (Estado, paTs, continente) e, 
portanto ele nao e muito adequado para o nTvel municipal. Para obter 
uma estimativa da area plantada com a cultura da cana-de-agucar na 
regiao de estudo, seria necessario que a escolha dos segmentos fosse 
aleatoria, o que na verdade nao ocorreu neste estudo cujo principal 
objetivo foi determinar o desempenho do procedimento para estimativa 
da proporgao da cultura nos segmentos teste. 

2.3 - METODOLOGIA PARA ESTIMATIVA DE PRODUTIVIDADE 

2.3.1 - MODELO AGROMETEOROLOGICO 


0 modelo agrometeorologico visa determinar a influencia 
da agua disponivel para a cultura sobre o seu rendimento. 

0 desenvolvimento do modelo baseou-se principalmente no 
relatorio numero 33 da FAO - Yield Response to Water, no qual esta 
descrita uma metodologia para calcular o renuimento maximo e 
rendimento real da cultura em fungao da disponibil idade de agua. 

A quantificagao do efeito do estresse hTdrico sobre o 
rendimento da cultura e obtido relacionando o decrescimo do rend.mento 
relativo (1 - RR/RM) e o deficit de evapotranspiragao relativa (1 - 

ETR/ETM), atraves do fator de resposta ao rendimento (Kr), por meio da 
seguinte equagao: 


(1 - RR/RM) 


Kr (1 - ETR/ETM), 


(2.3) 
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onde: 


RR = Rendimento Real , 

RM = Rendimento Maximo, 

ETR = Evapotranspiracao Real, 

ETM = Evapotranspiracao Maxima, 

Kr = Fator de resposta ao rendimento. 

Quando a quantidade d>i agua disponTvel for suficiente e 
suprir adequadamente a demands de agua pela cultura, entao ETR = ETM 
e, consequentemente, o rendimento maximo nao sera penal izado. No 
entanto, quando a demanda de agua pela cultura for maior do que a agua 
disponivel, entao ETR < ETM e, consequentemente, RR < RM. 0 fator Kr 
foi obtido empiricamente e se encontra tabelado em Doorembos e Kassam 
(1979). Para a cultura da cana-de-acucar o fator Kr e igual a 1,2 o 
que implica num deficit de rendimento relativo proporcionalmente, 
maior do que o deficit de evapotranspiracao relativa. 

0 rendimento maximo (RM) para a cultura da 
cana-de-acucar na regiao de Aracatuba e Guararapes foi calculado 
baseado no metodo desenvolvido por Kassam (Doorembos e Kassam, 1979). 
Os fatores climaticos que determinam o valor de RM sao temperatura e 
radiacao, alem do numero de meses do ciclo da cultura. E considerado 
que condicoes adversas tais como agua, nutrientes, pragas, doencas, 
ervas daninhas e outros nao sejam limitantes para o crescimento da 
cultura. 


No calculo da evapotranspiracao maxima (ETM) utiliza-se 
0 conceito da evapotranspiracao de referencia (ETO), a qual e 
relacionada com a ETM atraves do coeficiente da cultura (Kc) que 
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ajusta os valores de ETO para ETM de acordo com o estagio de 
desenvolvim^nto da cultura. Para o calculo de ETO emprega-se a equacao 
de Penman. 


A evapotranspiraqao real (ETR) e obtida a partir de urn 
Tndice de agua disponTvel no solo que utiliza dados de precipitacao, 
capacidade de retencao de agua do solo e a quantidade de agua 
remanescente no solo calculada por meio do balanqo hTdrico de 
Thornthwaite-Mather (1955). Quando este Tndice e maior ou igual a 1 
entao ETR = ETM, caso contrario utiliza-se uma tabela para obter os 
valores de ETR com as informacoes acima mencionadas. 

Finalmente, pode-se obter o valor de rendimento real 
(RR) resol vendo a equacao: 


RR = RM 


1-K^ 0 - 


ETR j 
ETM 


(2.4) 


Todos os calculos do model o sao feitos com medias 
mensais dos dados meteorologicos. 

2.3.2 - OBTENCAQ DE INDICES VEGETATIVOS PARA ESTIMAR RENDIMENTO DA 
CANA-DE~ACOCAR 


Para a estimativa da produtividade atraves de indices 
vegetativos obtidos a partir de imagens LANDSAT, utilizando o 
IMAGE-100, requer-se que seja investigada a correlacao destes indices 
com a produtividade da cultura, medida no campo, em diversas areas. 

Primeiramente, a partir de fotografias aereas coloridas 
falsa-cor, na escala de 1:30.000 e obtidas pelo INPE em junho de 1983, 
foram selecionados 42 talhoes de cana-de-aqucar e, a seguir, plotados 
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sobre as cartas topograficas na escala de 1:250.000. Munidos destes 
documentos cartograficos, foram coletados em campo sobre os talhoes 
selecionados dados referentes a sua area e producao, variedades, datas 
de plantio e colbaita da canaj discriminando para esta ultima 
infonnacao se a cana era de 19, 29 ou 39 corte. 

Em segundo lugar, a parti r de dados do LANDSAT, em fitas 
CCTs (fitas compatTveis com o computador), referentes as datas de 
passagem de 09/set. /82, 28/nov./82, 16/fev./83, 07/maio/83, 
10/jul./83, 26/jul./83, 27/ago. /83 e 28/set. /83, foram obtidas nos 
quatro cana is do MSS do LANDSAT informacoes espectrais de cada urn dos 
42 talhoes de cana selecionados, atraves do Analisador de Imagens 
Multiespectrais (IMAGE-100). 0 programs utilizado no sistema IMAGE-100 
para obtencao das informacoes espectrais foi o SINGLE-CELL. 

Devido a problemas de cobertura de nuvens que ocorreram 
nos dados referentes as passagens do LANDSAT de 28/nov./82 e 
07/maio/83, foi impossTvel obter informacoes espectrais de alguns 
talhoes de cana. Os dados referentes a passagem de 28/set. /83 foram 
abandonados em virtude dos talhoes de cana apresentarem ao sensor MSS 
uma forte influencia do solo em relacao a massa vegetal, ja que este 
perTodo corresponde, praticamente, ao fim da safra da cana-de-acucar. 
Alem destes problemas, nao foi possTvel a obtencao de inforrnacao de 
campo sobre o rendimento de todos os 42 talhoes de cana selecionados. 

Os 24 talnoes restantes eram compostos de cana planta e 
cana soca e que na maior parte haviam sido plantados em abril de 1980 
e abril de 1981, respectivamente. 

Para a analise da correlacao simples entre o rendimento 
(ton/ha) e os dados espectrais da cana em cada urn dos canais do 
LANDSAT, foram utilizadas somente as passagens de 16/fev./83, 
07/maio/83, 10/jul./83, 26/jul./83 e 27/ago. /83. A partir destes dados 
foram obtidos tambem os seguintes indices vegetativos: 
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D ; diferen^a de contagem digital entre o canal 5 e o canal 7, 

5-7 

% 

R : razao de contagem digital do canal 7 sobre o canal 5, 

7-5 

IVT : Tndice vegetative transformado, calculado atraves da se 
guinte equagao: 


IVT 


C6 - C5 
C6 + C5 


0,5 


(2.5) 


onde: 


C5 e C6 correspondem a respostas nos canais 5 e 6 do MSS uo 
LANDSAT, respectivamente, ■ 

IVT ; Tndice vegetative transformado, obtido atraves desta 

7 

equagao: 

IVT = 0^ (2.6) 

7 / C7 + C5 

Para testar a significancia dos coeficientes acima foi 
utilizado urn a = 0,05. 
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2.3.3 - PRODUTIVIDADE BASEADA NA AREA DE CANA CORTADA 


A principal razao que levou a realizaciao desta atividade 
esta fundamentada no fato de que e possTvel ter o controle diario do 
total de cana cortada que alimenta uma unidade industrial de acucar ou 
de alcool . No caso especTfico dos municTpios de Aracatuba e Guararapes 
esta condicao e mais favorecida, pois so existem 5 destilarias em 
operacao. 


Para a realizacao desta atividade foram utilizadas 
imagens LANDSAT do canal 5, na escala de 1:250. OuO, referentes as 
seguintes passagens: 10 e 26/jul./83, 27/ago. /83 e 28/set. /83. 
Utilizaram-se tambem aerofotograf ias coloridas falsa-cor, na escala de 
1:30.000, obtidas em junho de 1983, que cobreu: toda a area dos dois 
municipios (Figura 1.1). 

A obtencao dos dados de producao de cana foi feita 
atraves de trabalho de campo junto as destilarias Alcoazul, Aralco, 
Cruzalcool, Destivale e Unialco. 

Como metodo utilizado para a execucao desta atividade foi 
feita, inicialmente, a plotagem de todas as areas de cana sobre a 
passagem de 10/jul./83, utilizando fotografias aereas infravermelhas 
coloridas obtidas pelo INPE. Avaliou-se, tambem, a area de cana cortada 
sobre estas fotografias. 

Em seguida, sobre cada imagem do LANDSAT disponivel, 
fez-se a avaliacao de area do tema cana cortada atraves de 
interpretacao visual do canal 5 correspondente a regiao do vermelho do 
espectro eletromagnetico. 

Para a quant1fica<:ao do tema cana cortada utilizou-se 
grade milimetrada. 


- 31 - 


Os dados de tonelagem de cana inoTda foram obtidos junto 
as destilarias nas mesmas datas das passagens do LANDSAT e no periodo 
corresponden ao aerolevantamento. 

A razao entre a tonelagem de cana moTda e a area de 
cana cortada, ate uma determinada data, forneceria o rendimento medio 
da cultura obtido para as cinco destilarias desses municTpios. 


CAPITULO 3 


RESULTADOS E DISCUSSAO 


3.1 - DADOS OBTIDOS ATRAVES DE FOTOGRAFIAS AEREAS 

3.1.1 - distribuicao espacial das propriedades agricolas 


Nos mapas do atlas anexo pode-se observar a distribuicao 
espacial das propriedades agricolas levantadas dentro da area piloto. 

A Tabela 3.1 mostra o numero de propriedades agricolas 
que cultivam a cultura da cana-de-acucar, financiadas pelo Banco do 
Brasil S.A. nos dois municipios levantados. 

TABELA 3.1 

RESULTADOS OBTIDOS REFERENTES AO NOMERO DE PROPRIEDADES AGRICOLAS, QUE 
CULTIVAM A CANA-'DE-ACDCAR, FINANCIADAS PELO BANCO DO 
BRASIL S.A., DENTRO DA AREA PILOTO 


municTpio 

HESTILARIA 

NO DE PROPRIEDADES 

PERCENTAGEM 


ALCOAZUL 

22 

26 


ARALCO 

22 

26 

ARAQATUBA 

CRUZALCOOL 

11 

12 


DESTIVALE 

15 

18 

GUARARAPES 

UNIALCO 

15 

18 

TOTAL 

85 

100 


I'Am; i>i.AXK xcrr 
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A area de cada utna das 85 propriedades agricolas 
levantadas e ocupadas com cana-de-agucar foi obtida de duas maneiras: 
por informacao verbal de cada entrevistado e atraves de contagem de 
pontos de grade, sobre o mapa tematico, que contem os limites das 
areas, na escala de 1:250.000. 

A label a 3.2 compare os resultados dos valores de area 
e mostra as diferencas absolute e relative entre os dois processes de 
obtencao de area para cada uma das 85 propriedades. 

T ABEL A 3.2 

AREAS DAS PROPRIEDADES AGRiCOLAS PRODUTORAS DE CANA, OBTIDAS 
ATRAVES DE INFORMACAO VERBAL NO CAMPO E ATRAVCS DE CONTAGEM 
DE PONTOS SOBRE OS MAPAS TEMATICOS 1:250.000 


DCSTILAXIA 

COOICO DA 

PROPRIEDAOC ACRICOLA 

Xrca das propriedaocs agricolas (H«) 

DIFERENCAS 

informacao 

VERBAL (A) 

CONTAGEM DE 
PONTOS (B) 

8 . A 

® * A . 

B 


46 

1.189 .65 

1.206.25 

16.60 

1,40 


»1 

1.035.76 

1.100.00 

64.24 

6.20 



?79.52 

1.181.25 

201 ,73 

20,60 


52 

1,337.05 

1.293.75 

- 43.30 

- 3,20 


34 

2.262.70 

2.256.25 

- 6.45 

> 0.28 

•J 

26 

3.509.00 

3.606.25 

97,25 

2.77 

3 







43 

188.08 

381.25 

193.17 

102.70 

< 

SO 

629.20 

656.25 

27.05 

4.3C 

O 






u 

36 

713.32 

662.50 

* 50.82 

• 7.12 


37 

628.38 

631.25 

2.87 

0.45 

< 







3t 

526.45 

537.50 

11.05 

2,10 


41 

188.11 

187.50 

- 0.61 

- 0.32 


(continue) 
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Tabela 3.2 - Contlnua^ao 


ORIGIN.U PAGc !S 
OF POOR QUALITY 


DESTILARIA 

• 

cOdigo da 

PROPRIEOAOE AGRiCOLA 

AREA DAS PROPRIEDADES AGftICOLAS (Ha) 

DIFERENCAS 

informacao 
VERBAL (A) 

CONTAGEM DE 
PONTOS (B) 

8 - A 

8 - A .. 

B 


40 

— 

275.00 

— 

ARRENDAMENTO 


45-A 

1.189,00 

1.231,25 

42.25 

3.55 

•J 

42 

188,00 

206,25 

18,25 

9,70 

s 

47 

1.036,00 

1.262,50 

226,50 

17,94 








285 

2.541,00 

2.675,00 

134,00 

5,01 

< 






O 

66 

— 

787,50 

— 

— 

u 







38 

722.97 

675,00 

- 47,97 

• 6,60 







< 

35 

711,27 

725,00 

13.73 

1.93 


39 

— 

112,50 


ARRENDAMENTO 


30 

— 

3t3,75 

— 



18 

287,98 

281,25 

. 6,73 

- 2,39 


99 

704,22 

662.25 

- 41,97 

- 6,33 


22 

1.559.63 

1.650,00 

90.37 

5,47 


7 

1.381,82 

1.337,50 

- 44,32 

- 3,31 


1 

96,80 

100,00 

3,20 

3,20 

o 






u 

5 

609,84 

568,75 

- 41,09 

- 7,22 

u 

3 

1.04S.44 

1.068,75 

23.31 

2,18 

< 






ec 

2 

125,84 

131,25 

5.41 

4.12 

< 

20 

121,00 

187,50 

66,5 

35.47 


4 

774,40 

750,00 

- 24.40 

- 3,25 


6 

338,80 

281.25 

- 57,55 

- 20.46 


8 

208,12 

231,25 

23,13 

10,00 


68 

— 

47S.OO 

— 

ARRENDAMENTO 


(continua) 
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label a 3.2 - Conti nuagao 


OF . - 


CODIGO DA 

PROPRIEOAOE ACRlCOLA 



Area das proprieoaoes agrTcolas (Ha) 

D2FERENCAS 

ikformacao 

VERBAL (A) 

CONTAGEM DE 
PONTOS (B) 

■ 

° ‘ * * 100 
8 

1.887,60 

1.918,75 

31,15 

1.62 

486,42 

462,50 

- 23,92 

- 5,17 

246,84 

225,00 

- 21,84 

- 9,70 

171,82 

156,25 

- 15,57 

- 9.96 

242,00 

262,50 

20,50 

7,80 

847,00 

800,00 

- 47.00 

- 5,86 

1.210,00 

1.156.25 

- 53,75 

- 4.64 


4.850,0b 

27,91 

0,57 

262,10 

231,25 

- 30,85 

- 13,34 


3.872.00 

363.00 

363.00 

3.146.00 

1.573.00 
534,82 

1.210.00 
162,74 
266,20 
266,20 

2.037,64 


3.756,25 

387.50 
193,75 

3.662,50 
1.800.00 
906,25 
1.150 00 

137.50 
275,00 

312.50 
2.018,75 


- 3,08 
6,32 

- 87,35 
14,10 
12,60 
40,98 • 

- 5,21 

- 18,35 

3,20 

14,31 

- 0,93 












Tahela 3.2 - Conclusao 


OF POOR QUALii 


DCSTILARIA 


COOIGO DA 

PROPRUnAOC ACRICOLA 



*RtA DAS PROPRIEDADES AGRiCOLAS (Ha) 

oifekencas 

INFORMACAO 

CONTACEH OE 

n 

IB 

■ ^ V inn 

VERBAL (A) 

PONTOS IB) 

H 

B 

859,10 

906,25 

47.15 


5,20 

2.219.62 

2.518,75 

299,13 


11 .87 

380.30 

625,00 

244.70 


'39,15 

2.076.95 

2,056,25 

- 20.70 

- 

1,00 

290.16 

600,00 

309,84 


51.64 

256.58 

900,00 

643,42 


71,49 

292.75 

531.25 

238,80 


44,95 

751.31 

818.75 

67,44 


8.23 

A2.95 

43.75 

0,80 


1,82 

145.20 

187.S0 

42.30 


22,56 

683,13 

675,00 

- 8,13 

- 

1 ,^U 

583,20 

637,50 

54.30 


8,51 

3.400,10 

3.537.50 

137,40 


3,88 

314,60 

218,75 

- 95,85 

- 

43.81 

3,339,60 

3.412.50 

72,90 


2.13 

191,18 

212,50 

21,23 


10,03 

2.144,12 

2.225.00 

80,88 


3.63 

121,00 

150,00 

29.00 


19,33 

726,00 

781 ,25 

55,25 


7,07 

290,40 

306.25 

15,85 


5,17 

544,50 

525;00 

- 19,50 

- 

3.71 

3.339,60 

3.387,50 

47,90 


1.41 

2.178,00 

1.987,50 

-190,50 

- 

9,58 

431,58 

475,00 

- 6,58 

- 

1,38 

2.599,08 

2.681.25 

82,15 


3,06 

2.257,86 

2.350,00 

92,14 


:,92 
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Analisando a Tabela 3.2 pode-se cbservar que de 5 
propriedades agricolas (as de codigo 40, 66 ^ 39, 30 e 68) nao se 
conseguiu informacao' verbal sobre a area. Deste modo, considerando 
somente 80 propriedades agricolas, 16% dessas propriedades levantadas 
apresentaram uma discrepancia de valor de area da mais de 20%, entre os 
dois processos de obtencao. 

E possTvel que uma das maicres causas dessa percentagem 
de ocorrencia de discrepancia entre os valorem de area obtidos pel os 
;dois processos, mostrados na Tabela 3.2, tenha s'ia origem no erro de 
informacao verbal , a qucil foi obtida sem con provacao por documentos 
cartograficos ou juridicos. 

Na Tabela 3.3 e mostrado qi e houve uma discrepancia de 
4,51% entre o resultado obtido pela injrmacao verbal e pela contagem 
de pontos para as 85 propriedades agricolas. Esta tabela mostra ainda 
as discrepancias verificadas para cada destilaria. 

TABELA 3.3 


somatoriq da Area total, for destilar ia , das propriedades 

AGRlCOLAS PRODUTORAS DE CANA-DE-ACOCAR 


DESTILAR . 

Area total DASPROPR.AGRlCOLAS(Ha) 

DIFERENgA 

INFORNAgAO 
VERBAL (A) 

CONTAGEM DE 
PONTOS (B) 

B - A 

^xlOO 

ALC0A2UL 

19.575,46 

20.475,00 

899,54 

4,39 

ARALCO 

17.429,76 

17.312,25 

-117,51 

0,68 

CRUZALCOOL 

13.794,60 

14.600,00 

805,40 

5,52 

DESTIVALE 

9.964,03 

11.937,50 

1973,47 

16,53 

UNIALCO 

22.510,82 

22.887,50 

376,68 

1,65 

TOTAL 

83.274,67 

87.212,25 

3937,58 

4,51 
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3.1.2 - PRECISAQ CARTOGRAFICA OBTIDA NO TRACADO DOS LIMITES 
DAS PROPRIEDADES AGRlCOLAS 


Os mapas do atlas anexo, na escala de 1:50.000, e que 
contem os limites das propriedades agrTcolas, foram feitos por 
compilacao nao sendo realizados services de campo, tais como de 
topografia, de geodesia ou de reambulacao, para a determinacao de 
pontos de control e. 

UtiTizando o ZTS (Zoom Transfer Scope) procedeu-se a 
transferencia dos limites das propriedades dumarcadas nas fotografias 
aereas de 1:20.000 e 1:25.000 para as cartas topograficas na escala 
de 1:50.000. Por ocorrerem grandes diferencas entre a data das 
fotografias aereas e a da elaboracao das cartas, nem sempre foi 
possTvel urn ajuste perfeito; recorrendo, nestes casos, ao ajuste de 
pelo menos dois acidentes que pudessem definir uma orientacao para 
poder encaixar a figura desejada sobre a carta em 1:50.000. 

Tendo sido lancados os limites das propriedades 
agrTcolas sobre as cartas topograficas, impressas em "offset", na 
escala de 1:50.000, procedeu-se ao prepare do desenho cartografico 
sobre poliester de base estavel. 

A representacao grafica dos limites das propriedades 
agrTcolas foi uma fonte de erro, responsavel por uma parcela 
consideravel de imprecisao, a qual e tao mais acentuada quanto menor 
a area da propriedade agrTcola. Na Tabela 3.4 e mostrado, o erro que 
se comete, em funcao do tamanho da area, ao se utilizar para o 
desenho grafico uma pena de 0,1mm de diametro. 


TABELA 3.4 


ORIGINAL PAG.: .r 
OF POOR quality 


TABELA DEMONSTRATIVA DO ERRO DEVIDO A REPRESENTACAO GRAFICA AQ SE 
UTILIZAR PARA 0 DESENHO, PENA DE OJmm DE DIAMETRO 


Dimensoes da area em nm 

(Inn X Inm) 

(2nn X 2mm) 

(5mn X 5mm) 

(10mm X 10mm) 

(20mm X 20mm) 

Urea correspondente em 1:250.000 

6.25 ha 

25 ha 

156,25 ha 

625 ha 

2500 ha 

Rrea correspondente em 1: 50.000 

0,25 ha 

1 ha 

6,25 ha 

25 ha 

100 ha 

ERRO EM X 

i los 

i 5S 

i 2t 

CM 

V 

< 2X 


Como para o desenho em 1:50.000 neste trabalho foi 
usada caneta de 0,5mm de diametro, tem-se a seguinte Tabela 3.5 
comparativa; 


TABELA 3.5 

TABELA DEMONSTRATIVA DO ERRO DEVIDO A REPRESENTACAQ GRAFICA AO $E 
UTILIZAR PARA 0 DESENHO, PENA DE 0,5mm DE DIAI'IETRO 


Dimensoes da area em min 

(1mm X 1mm) 

(2mm X 2mm) 

(5mm X 5im) 

o 

X 

o 

(20im X 20na) 

Rrea correspondente em 1:250.000 

6,25 ha 

25 ha 

156,25 ha 

625 ha 

2500 ha 

Urea correspondente em 1; 50.000 

0,25 ha 

1 ha 

6,25 ha 

25 ha 

100 ha 

ERRO EM 2 

- 18% 

- 302 

t-t- 

o 

i 52 

LD 

V 


Comoo material foi todo trabalhado em 1 :50.000 e 
depois reduzido para a escala de 1:250.000, tem-se finalmente os 
seguintes valores da Tabela 3.6 




TABELA 3.6 


I 

TABELA DEMQNSTRATIVA DO ERRO DEV IDO A REPRESENTACAO GRAFICA AO 
UTULIZAR, PARA 0 DESENHQ, PENA DE 0,5mm DE DlAl^TRO NA ESCALA 
DE 1.50.000, REDUZIDA PARA 1:250.000 


Dimens5es da area em mm 
na escala de 1:50.000 

(1mm X lim) 

(2nri X 2nn) 

{5nm X 5mm) 

(10mm X lOmm) 

(20mm X 20mm) 

Ditnensoes da area em mm 
na escala de 1 :250.000 

(D,2nni x O.Zinn) 

(0,4mr: x 0,4trm) 

(1mm X 1mm) 

(2mm X Ztim) 

(4mm X 4mm) 

ERRO EM X 

indetenninado 

indetenninado 

t 10X 

i 5X 

- 2X 


Os dados da label a 3.6 demonstram que o tamanho mTnimo 
de uma propriedade, teoricamente viavel a ser mapeada, na escala de 
1:250.000, nao podera ser inferior a Imm^, equivalente a 6,25 ha, 
incorrendo-se era urn erro de 10% causado pela sua representagao 
graf ica. 


As imagens LANDSAT utulizadas neste trabalho 
apresentaram variagoes de escala muito reduzidas, sendo o raenor 
desvio medio encontrado igual a 0,33% e o maior igual a 0,85%, era 
relagao a carta topografica do IBGE. 
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3.2 - DADOS OBTIDOS ATRAVES DA INTERPRETACT^O DE IMAGENS DO LANDSAT 

« 

3.2.1 - ESTIMATIVA DA AREA DE CANA-DE>ACOCAR CULTIVADA POR 

PROPRIEDADE AGRlCOLA 

No atlas anexo i mostrada a distribuicao espacial da 
cultura canavieira obtida atraves da interpretagao visual das imagens 
LANDSAT. A este mapa e superposta uma transparencia com o resultado 
da interpretacao das fotografias aereas, com a final idade de 

comparacao dos dois resultados. 

A comparacao quantitativa destes resultados referentes 
a cana safra/83 e ao plantio/83 e mostrada na Tabela 3.7. 

TABELA 3.7 


RESULTADOS COMPARATIVOS DE AREAS OCUPADAS COM CANA-DE-ACPCAR, POR 
PROPRIEDADE AGRlCOLA, OBTIDOS A PARTIR DE DADOS DE AERONAVE E 

DADOS DO LANDSAT 


DESTILARIA 

cOdigo da 

— 
Area de cana (Ha) 

DIFERENCA (ha) 

PROPRIEDADE 

AERONAVE (A) 

LANDSAT (B) 

A- B 






46 


481,25 

- 118,75 

- 32,76 


51 

412,50 

112,50 


72,73 

=> 

M 

<I 

33 

18,75 

143,75 


- 666,67 

O 

o 

52 




31,25 

_J 

eC 

34 

837,50 



58,21 


26 

1.287,50 

737,50 


42,72 


(continua) 
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Tabela 3.7 - Conti.nua^ao 


ESTILARIA 


CODIGO DA 

Area de cana (Ha) 

PROPRIEDADE 

AERONAVE(A) 

LANDSAT(B) 

97 

162,50 

125,00 

65 

68,75 

75,00 

69 

131,25 

112,50 

27 

650,00 

431,25 

98 

187,50 

143,75 

19 

87,50 

87,50 


1.168,75 
5 


X 100 


23,08 

9,09 

14,29 

33,65 

23,33 

0,00 


45,99 

34,78 

50,00 

20,63 

29,58 

9,13 

100,00 


(continua) 
























DESTILARIA 

COD I GO DA 

AREA DE CANA (Ha) 

DIFERENCA (ha) 

PROPRIEDADE 

AERONAVE (A) 

LANDSAT (B) 

B 



316 

262,50 

331,25 

- 68,75 

- 25,19 


317 

843,25 

981,25 

- 138,00 

- 16,37 


305 

62,50 

0,00 

62,50 

100,00 


323 

362,50 

181,25 

181,25 

50,00 


320 

131,25 

50,00 

81,25 

61,90 


318 

656,25 

818,50 

- 162,25 

24,72 

o 






o 

319 

181,25 

43,75 

137,50 

75,86 

«£ 

310 

62,50 

50,00 

12,50 

20,00 


306 

81,25 

0,00 

81,25 

100,00 

Z 






ZD 

308 

106,25 

0,00 

106,25 

100,00 


321 

225,00 

125,00 

100,00 

44,44 


314 

112,50 

175,00 

- 62,50 

- 55,56 


324 

100,00 

31,25 

68,75 

68,75 


315 

918,75 

756,25 

162,50 

17,69 


325 


281,25 


0,00 


281,25 


100,00 












Com a final idade de apresentar a Tabela 3.7 de forma 
resumida, foi elaborada a Tabela 3.8 na qual os resultados 

comparat^ os de areas ocupadas com cana-de-acucar por propriedade 
agrTcola, obtidos a partir de dados de aeronave e dados do LANDSAT, 
sao apresentados por destilaria. 

TAB ELA 3.8 

RESULTADOS COMPARATiVOS DE AREAS OCUPADAS COM CANA-DE-ACOCAR POR 
PROPRIEDADE AGRTCOLA, OBTIDOS POR DESTILARIA A PARTP DE DADOS 

DE AERONAVE E DO LANDSAT 



NOMERO DE 

AREA DE CANA (Ha) 

DIFERENCA (Ha) 

DESTILARIA 

PROPRIEDADES 

AERONAVE (A) 

LANDSAT (B) 

A - 8 

xlOO 

ALCOAZUL 

22 

7.081 ,00 

4.806,25 

2.274,75 

32,12 

ARALCO 

22 

4.625,00 

3.387,00 

1.238,00 

26,77 

cruzalcool 

11 

3.576,75 

2.680,00 

896,00 

25,05 

DESTIVALE 

15 

4.968,75 

4.162,50 

806,25 

16,23 

UNIALCO 

15 

4.387,00 

3.543,50 

843,50 

19,23 

TOTAL 

85 

24.638,50 

18.580,00 

6.058,50 

24,59 
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Na analise da relacao entre os dados de area obtidos 
por imagens LANDSAT e por aeronave observa-se que, de um modo geral, 
as imagens tendem a subestimar os valores de area. Este fato pode ser 
verificado atraves do grafico da Figura 3.1. 

A razao entre a estimativa de area de cana-de-acucar 
obtida com imagens LANDSAT e aquela obtida com dados de aeronave foi 
de 75%. 


Com base nos dados da Tabela 3.7 foi obtida uma 
equacao de regressao que expressa a correlacao linear entre os 
valores de area de cana, por aeronave e por LANDSAT: 

? = 1,01X + 68,42. 

0 grafico da Figura 3.1 mostra a relaqao dos dados 
analisados e a respectiva reta de regressao. 

0 coeficiente de correlacao (r = 0,88) foi 
significative ao nivel de 1% de probabil idade. 

A interpretacao das imagens LANDSAT, utilizando varias 
passagens do satelite, ficou prejudicada devido a cobertura de nuvens 
e a ausencia de informacao causada por problemas no sensor MSS do 
LANDSAT- 3, conforme mostra o mapa da Figura 3.2. 

Em decorrencia de tais fatos fez-se uma selecao das 
propriedades que nao foram afetadas pela falta de dados e gerou-se 
uma nova analise de regressao, utilizando apenas os resultados de 36 
propriedades agrTcolas, cuja equaqao e; 


Y = 0,90X + 95,00. 


AREAS OBTIOAS COM DADOS DE AERONAVE (Hd 







ORIGINAL PAG£ iS 
OF POOR QUALITY 



Fig. 3.2 - Caracterizagao de areas com problemas de cobertura de 
nuvens e perda de informafio em imagens do LANDSAT. 
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A Figura 3.3 mostra a relacao dos dados selecionados 
pelo criterio acitna mencionado e a reta de regressao. 0 teste t 
mostrou que o coeficiente de correlacao (0,91) foi significative ao 
nTvel de 1% de probabil idade. 

A diferenqa entre os resultados obtidos da 
interpretacao dos dados de aeronave (A) e dos dados do LANDSAT (B), 
apresentados na Tabela 3.7, ten valor positivo quando a area de cana 
interpretada nas imagens LANDSAT foi menor do que aquela obtida pelos 
dados de aeronave. 

As diferengas positives entre A e B indicam que as 
areas de cana obtidas da interpretacao das imagens LANDSAT foram 
subestimadas. 


A fitn de avaliar os erros de omissao e os erros de 
inclusao sobrepds-se a interpretacao de cana das imagens do LANDSAT a 
interpretacao de cana das fotografias aereas, obtendo-se um mapa 
mostrado no Apindice A que contem a distribuicao das areas de 
cana-de-acucar interpretadas corretamente e os erros de omissao e 
inclusao. 

A partir desse mapa foi possTvel calcular, as areas 
interpretadas corretamente, as areas de cana nao interpretadas e as 
areas interpretadas incorretamente como cana nas imagens do LANDSAT. 
Na Tabela 3.9 sao mostrados, por propriedade, os valores das 
percentagens de cana interpretadas corretamente e as percentagens dos 
erros de inclusao e omissao, conforme as seguintes formulas; 


AREAS OSTIOAS COM DADOS DC ACRONAVC (N«) 



Fig. 3.3 - Regressao entre os valores da area de cana, 
obtidos com dados do LANDSAT e dados de 
aeronave, referentes as 36 propriedades 
agrTcolas que nao foram afetadas pela falta 
de dados do LANDSAT. 
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Pa = X TOO; El = a TOO; EO - 100 - Pa, 

ACR ACR 

« 

onde: 

Pa - Cana Interpretada Corretamente (%), 

El - Erro de Inclusao (%), 

EO - Erro de Om’ssao (^), 

ACC - Area de Cana Interpretada Corretamente (ha), 

ACR - Area Real de Cana (ha), 

ACI - Area de Cana Interpretada Incorretamente (ha). 

TABELA 3.9 

RESULTADO DA INTERPRETAQAO DE CANA A PARTI R DE IMAGENS DO 
LANDSAT E ERROS COMETIDOS 


DESTILARIA 

cOdigo da 

PROPRIEDADE 

CANA 

INTERPRETADA 

CORRETAMENTE 

{%) 

ERROS (%) 

INCLUSAO 

OMISSAO 


46 

100,00 

25,97 

0,00 

=3 

M 

51 

33,87 

0,00 

66,13 

< 

O 

52 

76,12 

0,00 

23,88 

o 
1 

34 

56,48 

4,70 

43,52 

c 

26 

63,10 

5,40 

36,90 


(continua) 












Tabela 3.9 - Continuagao 


DESTILARIA 

COD I GO DA 

CANA 

INTERPRETADA 

ERROS (%) 

PROPRIEDADE 

CORRETAMENTE 

C%) 

INCLUSAO 

OMISSAO 


5 

48,14 

0,00 

51,86 


3 

68,09 

13,51 

31,91 


2 

100,00 

0,00 

0,00 


20 

0,00 

0,00 

100,00 


4 

79,55 

7,89 

20,45 


6 

14,71 

0,00 

85,29 

o 

8 

100,00 

50,00 

0,00 

c_> 

68 

0,00 

0,00 

100,00 


9 

23,25 

33,30 

76,75 

< 





DC 

73 

33,30 

0,00 

66,60 


77 

0,00 

0,00 

100,00 


97 

83,30 

11,80 

16,70 


65 

100,00 

37,50 

0,00 


69 

57,14 

46,60 

42,86 


27 

55,17 

2,00 

44,83 


98 

0,00 

76,33 

100,00 


19 

100,00 

0,00 

0,00 

t 

10 

37,68 

38,82 

62,32 

o 





c_> 

I 

13 

78,57 

0,00 

21,43 

< 





ZD 

Ct 

14 

80,00 

0,00 

20,00 

o 






(continua) 
















COD I GO DA 

1 CANA 

INTERPRETADA 

! 

ERROS (%) 

PROPRIEDADE 

CORRETAMENTE 

m 

INCLUSAO 

OMISSAO 

317 

i 

73,70 

45,70 

26,30 

305 

0,00 

0,00 

100,00 

323 

52,00 

38,10 

48,00 

320 

50,00 

0,00 

50,00 

318 

91,90 

38,30 

8,10 

319 

56,50 

7,10 

43,50 

306 

0,00 

0,00 

100,00 

308 

0,00 

0,00 

100,00 

321 

17,60 

57,10 

82,40 

314 

78,60 

38,90 

21,40 

324 

54,50 

0,00 

45,50 

315 

63,20 

36,30 

36,80 

325 

0,00 

0,00 

100,00 
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Com 0 objetivo de apresentar os dados da Tabela 3.9 de 
forma resumida, foi elaborada a Tabela 3.10 na qual os dados de 
percentagem de cana interpretada corretamente e as percentagens de 
erros de inclusao e omissao sao apresentados por destilaria. 

TABELA 3.10 

VALORES DE CANA INTERPRETADA CORRETAMENTE E ERROS DE INCLUSAO 

E OMISSAO POR DESTILARIA 


DESTILARIA 

NOMERO DE 
PROPRIEDADE 

CANA INTERPRETADA 
CORRETAMENTE 
(%) 

ERROS (%) 

INCLUSAO 

OMISSAO 

ALCOAZUL 

21 

70,47 

9,91 

29,53 

ARALCO 

22 

.57,33 

21,62 

42,67 

CRUZALCOOL 

IQ 

52,06 

26,03 

47,94 

DESTIVALE 

15 

68,22 

19,58 

31,78 

UNIALCO 

13 

58,38 

36,73 

41.62 . 

TOTAL 

81 

63,12 

20,55 

36,88 


Analisando os resultadjs da Tabela 3.10 verifica-se que 
houve melhor desempenho He interpretacao da cana-de-acucar da 
Jestilaria Alcoazul (70,47%) e menor ' desempe*'ho da destilaria 
Cr-uzalcool (52,06%). A percentagem de cana interpretada (orretamente 
na area abrangida pelas cinco destilarias foi de 63,12%. 
































No teste para determinar a exatidao do mapeamento, 
considerando as classes cana e nao cana, constatou-se que dos 93 

I 

pontos sorteados para verificacao, 26 foram interpretados 
incorretamente. Isto indica que o mapa produzido nao atendeu os 
criterios preestabelecidos conforme descritos na metodologia (secao 
2. 2. 3. 5). I-lantido o mesmo risco de 5% do consumidor e produtor, 
pode-se afirmar que o mapa produzido teve uma exatidao minima de 
63,50%. 


Alguns fatos contribuTram para a ocorrencia dessa baixa 
exatidao de mapeamento tais como: nas imagens LANDSAT, a cultura da 
cana-de-agucar foi identificada atraves de uma chave de interpretagao 
descrita na segao 2. 2. 3. 2, que considera as caracterTsticas 
espectrais e temporais do alvo. Esta chave de interpretagao procura 
identificar a cultura de cana dentro dos seus padroes normais de 
desenvolvimento e estabelece criterios rigidos na interpretagao. 
Assim sendo, a identificagao da cultura ficou prejudicada quando esta 
nao se apresentou com a tonal idade prevista na chave de 
interpretagao. Esta condigao atTpica de tonal idade da cultura nas 
imagens pode ser atribuida a fatores diversos, tais como: tombamento 
da cana, incidencia de ervas daninhas, cana bisada, cana com 
florescimento, mistura va-ietal, erosao, tratos culturais e manejo. A 
estes fctores sao somcdos os efeitos da excessiva precipitagao que 
ocorreu no primeiro semestre de 1983 que prejudicou a caracterizagao 
de algumas areas carjvieiras em termos da resposta espectral. Este 
fato tambem impediu que as areas de cana pudessem ser melhor 
diferenciadas das areas de pastagem, que constituem uma das 
principals fontes de erro na interpretagao. 

A escala de trabalho das imagens do LANDSAT constitui 
tambem uma fonte de diverging i a entre os dados gerados neste si sterna 
e aqueles obtidos por fotos aereas. 
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3.2.2 - ABORDAGEH ALTERNATIVA PARA IDENTIFICACAO DE CANA, AO NTVEL DE 
PROPRIEDADE AGRtCQLA. ATRAVES DE IMAGENS DO LANDSAT 


Em adicao a metodologia que visa substituir, 

totalmente, o processo de fiscalizagao, uma metodologia alternativa 
foi desenvolvida com a finalidade de complementar a fiscalizaqao 
convencional (visita do fiscal ao campo), fornecendo informaqoes 
objetivas de avaliaqao de area e acompanhamento do desenvol vimento de 
cultures agricolas, em propriedades que obtiveram f inanciamento. 

Os requisites necessarios para o alcance deste 
objetivj foram os seguintes: 

1) cadastre das propriedades agricolas, numa base cartografica 
na escala de 1:250.000; 

2) informaqoes referentes a area plantada com a culture da 
cana-de-agucar, destacando a area financiada e calendario de 
eventos, principalmente, data de plantio; 

3) croquis do uso da terra na propriedade agricola e 
estimative, em percentagens, dos diversos uses da terra 
(pastagem, mate, cultures temporaries e permanentes, etc). 

Em seguida, adotou-se uma abordagem metodologica que 

consiste em: 


1) obter as informacoes, junto ao Banco do Brasil, de areas e 
dates de plantio de cana nas propriedades agricolas 
investigadas; 

2) selecionar imagens do LANDSAT coincidentes com a epoca de 
plantio e delimiter areas com padroes de solo preparado; 
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3) verificar passagens subsequentes do LANDSAT e tentar 
identificar areas realmente plantadas com a cultura da cana 
sobra a^ areas delimitadas no Hem ?; 

4) delimitar, concomitantemente aos itens 2 e 3, todas as areas 
que potencialmente podem ser interpretadas como cana, alem 
de delimitar novas areas preparadas para plantio. 

As atividades referentes aos itens 2, 3 e 4 sao 

interativas, ou seja: a identificacao, em imagens do LANDSAT, de 

areas plantadas com cana, como areas delimitadas com alta 
probabilidade de ser cana foi feita em toda a sequencia de imagens 
utilizadas. 

Entretanto, a ausencia de imagens de boa qualidade 
referentes a perTodos considerados ideais para a identificacao da 
cultura canavieira interferiu nos resultados da interpretacao. Este 
fato, portanto, impossibilitou que algumas areas inicialmente 
identificadas como solo preparado ou mesmo com algumas 
caracterTsticas de talhoes de cana, nao pudessem ser confirmadas como 
area de cultura canavieira. A confirmacao destas areas deve ser feita 
posteriormente atraves de trabalho de campo executado pelo fiscal. 
Deve-se ressaltar que na presente metodologia o trabalho de campo foi 
substituTdo pela utilizacao de fotografias aereas. 

As areas da cultura canavieira que nao foram 
delimitadas como possiveis areas de cana e nem identificadas como 
cana, atraves das imagens, const ituem o erro de omissao da 
interpretacao. Por outro lado, as areas identificadas como cana nas 
imagens, mas que na realidade nao sao cana, constituem o erro de 
inclusao. A avaliacao destes dois tipos de erro devera ser feita 
utilizando o mapa da area canavieira obtido atraves de fotos aereas. 
As fotos aereas deverao ser utilizadas para tentar explicar a causa 
desses erros. 
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A necessidade de utilizar esta metodologia surgiu pelas 
seguintes razoes; os municTpios de Aragatuba e Guararapes nao tem, 
ate 0 mome’nto, tradicao de cultivo de cana-de-acucar, sendo uma 
regiao caracterizada pela exploracao da pecuaria. A implantacao de 
destilarias de alcool se iniciou recentemente (1980). Atraves de 
exaustivos trabalhos de campo na area piloto foram observados muitos 
problemas com a cultura, tais como: cana bisada, tombada e invadida 
por ervuZ daninhas. 

A seguir sao apresentados os resultados obtidos 
utilizando esta metodologia. 

A Tabela 3.11 resume os resultados da estimative de 
area, percentagem de area de cana interpreted? corretamente e 
percentagens dos erros de omissao e inclusao cometidos por 
destilaria ao adotar a abordagem descrita acima. 

Observa-se pelos resultados da Tabela 3.11 que a 
interpretacao de imagens do LANDSAT, com suporte em informacoes 
obtidas de campo se constitui numa abordagem altamente eficiente na 
avaliagao da area da cultura da cana-de-agucar, melhorando, 
substancialmente, a perceritagem de cana interpretada corretamente. 

Nota-se que houve melhor desempenho de interpretagao da 
cana-de-agucar na destilaria Destivale (97,95%) e o menor desempenho 
na destilaria Cruzalcool (88,95%). A percentagem de cana interpretada 
corretamente na area abrangid?. pei;s cinco destilarias foi de 91,77%. 

0 erro de inclusao (6,60% - areas interpretadas como 
cana, mas que, na realidade, sao outros alvos) teve como causa 
principal a identificagao de pastagens e de algumas cultures anuais, 
principalmente tomate e feijao, como cana-de-agucar. Ja o erro de 
omissao (8,23% - areas de cans nao interpretadas em imagens do 
LANDSAT), e explicado devido, principalmente, a ocorrencia de talhoes 
com cana bisada, tombada e com florescimento. 
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TABELA 3.11 

ESTIMATIVA DE AREA DE CANA, PERCENTAGEM DE CANA INTERPRETADA 
CORRETAMENTE E PERCENTAGENS DE ERROS DE OMISSAO 
E INCLUSAO FOR DESTILARIA 


DESTILARIA 

nOmero de 

PROPRIEDADES 

Area de cana interpretada 

CORRETAMENTE ATRAVES DE 
IMAGENS 00 LANDSAT SOM SU 
PORTE DE INFORMAgOES AUXT 
LIARES (Ha). 

Area real 

DE 

CANA 

(Ha) 

Area de cana 

INTERPRETADA 

CORRETAMENTE 

(i) 

ERROS (S) 

OMISSAO 

INCLUSAO 

ALCOALZUL 

20 

5.762,50 

6.406,25 

89,95 

10,05 

2,44 

ARALCO 

22 

3.931,25 

4.300,00 

91,41 

8,58 

6,25 

CRUZALCOOL 

10 

3.118,75 

3.506,25 

88,95 

11,05 

4,81 

DESTIVALE 

13 

3.875,00 

3.956,25 

97,95 

2,05 

6,00 

UNIALCO 

15 

2.756,25 

3.018,75 

91.30 

8,70 

18,84 

TOTAL 

84 

19.443,75 

21.187,50 

91,77 

8,23 

6,60 


Na analise das fotografias aereas observou-se, tambem, 
que os talhoes de cana localizados nas proximidades de cursos d'agua 
apresentaram tonalidades atTpicas e, portanto, nao interpretadas como 
cana-de-agucar. 


0 fato de que as imagens LANDSAT nem sempre puderam ser 
adquiridas durante os perTodos de melhor caracteri zagao da cultura.em 
toda a area piloto, deve ter contribuTdo, tambem de forma sign^ 
ficante, para a compos igao do erro de omissao. 

0 mapa do Apendice A mostra a area de cana interpretada 
corretamente e as areas com erros de omissao e inclusao cometidos na 
interpretagao das imagens do LANDSAT com o suporte de informagoes 
auxi Hares de campo. 

ORIGINAL PmZ IS 
OF POOR QUALin 



ORIGINAL P. c.- • ' - 

OF POuP C 

A Tabela 3.12 apresenta, por propriedade, a proporgao 
da area ocupada com cana-de-agucar, ohtida utilizando a presente 
abordagem. Como pode ser observado nesta tabela, o uso de urn grande 
numero de informagoes quer ao nivel orbital, quer ao nTvel de campo, 
nao foi suficiente para caracterizar totalmente a cultura da cana, em 
algumas propriedades. 


TABELA 3.12 

PROpORgAo DA Area ocupada com cana-de-acOcar, mo nivll de 

PROPRIEDADE, 03TIDA ATRAVES DE IMAGENS uO LANDSAT E COM 
TRABALHO SIMULADO DE CAMPO, PERCENTAGEM DE CANA 
I NTERPRETADA CORRETAMENTE E DETERMINACAO DOS 
ERROS DE OMISSAO E INCLUSAO COMETIDOS 


DESTILARIA 

cOdigo da 

X DE Area de 
cana obtida 

ATRAVES DE 

s DE Area de 

CANA CONE I R 
MADA POR TRA 
BALHO SIMULA 
DO DE CAMPO 

Area de cana 
INTERPRETADA 

ERROS (X) 

PROPRIEDADE 

FISCAlIZAChO 
COM IMAGENS 
00 LANDSAT 

CORRETAMENTE 

(2) 

OMISSAO 

INCLUSAO 


46 

100,00 

0,00 

100,00 

0,00 

0,00 


51 

48.33 

36,67 

85,00 

15,00 

0,00 


5Z 

70,97 

20,97 

91.94 

8,06 

0,00 


34 

16,54 

71,43 

87,97 

12,03 

0,00 


26 

30,29 

56,25 

86,54 

13.46 

2,40 

-J 

43 

32,43 

51,35 

83,78 

16,22 

0,00 


36 

0,00 

66,67 

66.66 

33,33 

0,00 

M 

37 

38,89 

30,56 

69,44 

30,56 

0,00 

O 

31 

79,41 

8,82 

88,24 

11,76 

41,18 

o 

41 

36,84 

52,63 

89,47 

10,53 

0.00 

-J 

40 

50,00 

50,00 

100,00 

0,00 

0,00 

«c 

4 5- A 

94,32 

5,68 

100,00 

0,00 

0,00 


42 

70,59 

29.41 

100,00 

0,00 

0,00 


47 

90,48 

0,00 

90,48 

9,52 

0,00 


285 

79,41 

20,59 

100,00 

0,00 

17,65 


66 

100,00 

0,00 

100,00 

0,00 

0,00 


(continua) 
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Tahela 3.12 - Contlnua^ao 


OF PC 


DESTILARIA 

CODIGO DA 

X DE AREA DE 
CANA OBTIDA 
ATRAVES DE 

X DE area OE 
CANA CONFIR 

:lada por trS 

BALHO SIMULA 
DO DE CAMPa 

Area de cana 
interpretada 
CORRETAMENTE 

w 

ERROS (X) 

PROPRIEDADE 

FISCAL I ZACAO 
COM I MAGE NS 
DO LANOSAT 

OMISSAO 

INCLUSAO 


38 

66.67 

26,19 

92,86 

7.14 

0,00 

53 

ISI 

< 

O 

35 

43,75 

25,00 

68,75 

31.25 

0,00 

< 

39 

100,00 

0,00 

100,00 

0,00 

0,00 


30 

100,00 

0.00 

100.00 

0.00 

0,00 


18 

57,69 

0.00 

57,69 

42,31 

0,00 


99 

70,00 

25.00 

95,00 

5,00 

0,00 


22 

100,00 

0,00 

100,00 

0,00 

46,15 


7 

66,67 

17,46 

84,13 

15,87 

3,17 


1 

92,31 

0,00 

92,31 

7.69 

7,69 


5 

84,78 

0.00 

84,78 

15,22 

0,00 

o 

3 

70,42 

28,17 

98,59 

1.41 

2.82 

u 

2 

100,00 

0,00 

100,00 

0,00 

0,00 

«c 

20 

0,00 

0,00 

0,00 

100,00 

0,00 

ec 

4 

91,18 

0,00 

91,18 

8,82 

0,00 

< 

6 

81,58 

7,89 

89,47 

10,53 

5,26 


8 

100,00 

0,00 

100,00 

0,00 

0,00 


68 

0,00 

100,00 

100,00 

0,00 

0,00 


9 

88,89 

11,11 

100,00 

0,00 

20,37 


73 

100,00 

0,00 

100,00 

0,00 

0,00 


77 

100,00 

0,00 

100,00 

0,00 

42,86 


(continua) 
















Tabela 3.12 - Contlnuagao 


OESTILARIA 

COOIGO OA 

% DE AREA DE 
CANA OBTIOA 
ATRAVES DE 

X OE AREA DE 
CAHA CONE I R 
MADA POR TRA 

Area de cana 
interpretada 
corretamente 

ERROS (X) 


PROPRIEDADE 

h ISCAl IZACAu 
COM IMAGENS 
DO LANDSAT 

BALHO SIMULA 
DO DE CAMPO 

(2) 

OMISSAO 

INCLUSAO 


97 

40,00 

40,00 

80,00 

20,00 

0,00 


65 

50,00 

25,00 

75,00 

25,00 

50,00 

o 

69 

64,71 

23,53 

88,24 

11,76 

0,00 

< 

27 

59,63 

39,45 

99,08 

0,92 

5,50 

«c 

98 

92,86 

0,00 

92,86 

7.14 

0,00 


19 

76,92 

0,00 

74,92 

23,08 

0,00 


10 

78,16 

16,67 

94,83 

5.17 

2,87 


13 

100,00 

0,00 

100,00 

0,00 

27,78 

o 

14 

100,00 

0,00 

100,00 

0,00 

0,00 

o 

<_> 

16 

70,23 

9,92 

80.15 

19,85 

12,98 


17 

64,20 

16,05 

80,25 

19,75 

0,00 

c£ 

67 

62,07 

20,69 

82,76 

17,24 

0,00 

Nl 

103 

100.00 

0,00 

100,00 

0,00 

0,00 


105 

85,71 

0,00 

85,71 

14,29 

0,00 

O 

107 

100,00 

0,00 

100,00 

0,00 

0,00 


80 

73,40 

22,34 

95,74 

4,26 

0,00 


93.48 
96, 
100,1 


(continua) 
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Tabela 3.12 - Conclusao 


OESTILARIA 

COO I GO DA 

% DE area de 
CANA OBTIDA 
ATRAVES DE 
FISCAL IZACAO 
COM IMAGENS 
DO LANDSAT 

% DE Area oe 
CANA C0NE!R 

Area de cana 

INTERPRCTADA 

[ ERROS (Z) 

PROPRIEOAOE 

MADA POR TRA 
CALHO jiMUlA 
DO DE CAMPC 

CORRETAMENTE 

w 

OMISSSO 

INCLUSAO 


321 

0,00 

75,68 

75. G8 

24,32 

13,51 

o 

o 

314 

0,00 

88,89 

88,89 

11,11 

66,67 

—1 

< 

324 

0,00 

100,00 

100,00 

0,00 

0,00 


315 

58,68 

41,32 

100,00 

0,00 

14,88 

=3 

325 

100,00 

0,00 



100,00 

0,00 

35,90 


3.2.3 - TAMANHO MlNIMO DE AREA DA CULTURA DA CANA-DE-ACHCAR 
DETECTAVEL EM IMAGENS DO LANDSAT 


Quanto a area minima da cultura da cana-de-acucar 
possTvel de ser detectada em imagens do LANDSAT, observou-se que 
areas menores que 25 ha tornam-se dificeis de ser mapeadas na escala 
de 1:250.000. De acordo com Bahr (1984), o erro relativo que se 
comete na avaliacao de alvos em imagens do LANDSAT e de lOX para 
areas de 10 ha, o que vai depender do contraste apresentado, ou seja, 
quanto maior o contraste, maior a chance de deteccao precisa do alvo, 
dada a melhor definiqao do seu contorno. 

3.2.4 - AVALIAQAO DO USO DA TERRA DA AREA PILOTO 

Foi elaborado urn mapa de uso da terra para a area dos 
municTpios de A^aqatuba e Guararapes. 



Neste mapeamento foram utilizadas imagens do 
MSS/LANDSAT, referentes aos canais 5 e 7, da passagem de 28 de 
setembro de 1983, na escala de 1:250.000. 

Para a interpretacao das imagens fez-se uso de uma 
grade quadn'culada (2,5ram x 2,5mm), representando, no terreno, area 
de aproximadamente 39,0 ha, sendo esta abordagem uma modificaoao da 
proposigao de Genderen et alii (1978). 

A interpretacao foi realizada atraves da identificacao 
de elementos da imagem, tais como: tonal idade, forma, homogeneidade, 
etc. 


Em cada quadrTcula analisada foi considerado o tipo de 
uso predoiTiinante, isto e, aquele que ocupasse area superior a 50% da 
referida quadricula. 

Na identificacao das classes de uso da terra 
utilizou-se o procedimento proposto por Anderson et alii (1979), que 
se refere ao sistema de uso da terra, utilizar.do dados de 
sensor iamento remote ao nivel orbital . 

Neste trabalho foram definidas as seguintes classes de 
uso da terra: pastagem, areas agrTcolas, vegetacao arborea e areas 
urbanas. 


Os resultados referentes as classes identificadas e 
mapeadas encontram-se na Tabela 3.12. 

% 
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TABFLA 3.13 

USO DA TERRA DA AREA PHOTO OBTIDO ATRAVES DE IMAGENS DO LANDSAT 


TEMAS 

PASTAGEM 

Areas 

vegetacao 

Areas 

TOTAL 

VAL0RES\^^ 


AGRI COLAS 

arborea 

URBANAS 


Area (kmM 

2.825,00 

450 ,39 

42,96 

28,12 

3.347,26 

PERCENTAGEM 

84,4 

J 

13,5 

1.3 

0.8 

100,0 


Os dados da Tabela 3.13, representam apenas os valores 
obtidos para as classes ja mencionadas, nao incluindo areas ocupadas 
pelos corpos d'agua, areas de varzeas e pequenos nucleos urbanos. 

Atraves desta tabela, pode-se observar a nTtida 
predominancia da classe pastagem que representa 84,4% da area de 
estiido, denionstrando ainda a sua vocacao para pecuaria extensiva. 

0 mapa do uso da terra obtido encontra-se no 

Apendice B. 

3.3 - CLASSIFICACAO AUTOMATICA 

3.3.1 - RESULTADOS DA CLASSIFICACAO AUTOMATICA 

0 algoritmo de classificacao automatica, denominado 
Procedimento INPE, forneceu como resultados os seguintes produtos: 

a) urn mapa (printout) com a distribuicao espacial da 

classificacao automatica da cana-de-acucar , nos cinco 
segmentos testes; 
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h) a propor?ao da cultura da cana nos segmentos testes, obtida a 
parti’r da classificagao automatica; e, 

c) a proporgao da cultura nos segmentos testes, corrigida de even 
tuais tendencias ocorndas durante o processo de classificagao 
automatica. 

0 mapa da classificagao foi utilizado para a tomada de 
decisoes sobre os resultados da classificagao automatica de cana de ca^ 
da segmento, repetindo-se ou nao todos os passos do Froaedimento INPE 
se os resultados avail ados no referido mapa eram aceitos ou rejeita 
dos. Utilizou-se, tambem, o mapa de classificagao para determinar a 
percentagem de acerto, atraves da classificagao automatica de cana-de 
-agucar. 

Na Tabela 3.14 sao mostrados os resultados em termos de 
percentagem de acerto e proporgao da cultura, obtidos a partir do Pro 
cedimento INPE^ alem da proporgao da cultura obtida pela verdade te£ 
restre para fins comparatives. 


TABELA 3.14 


RESULTADOS DO PROCEDIMENTO-INPE 


SEGMENTO 

TESTE 

mi 

PROPORCSO DA 
CULTURA OBTI 
DA PELO CLAI 
SIFICADOR “ 
(%) 

PROPORCAO DA 
CULTURA COR 
RIGIDA (A) “ 

(%) 

PROPORCAO DA 
CULTURA OBTI 
DA ATRAVES 
DE VERDADE 
TERRESTRE (B) 

DIFERENCA 

RELATIVA 

(° ■ ^) X 100 

B 

1 

65,0 

13,6 

12,9 

11,9 

- 4 

2 

64,1 

13,8 

11.4 

11,8 

3 

3 

62,8 

20,4 

16,2 

15,8 

- 3 

4 

61,7 

17,8 

15,8 

18,6 

15 

5 

45,5 

10,0 

8.2 

7.9 

- 4 

media 

59,8 

15,1 

12,9 

13,2 

2 


I 


k 












I 

I 

1 

1 

I 

I 

I 

1 

I 

I 

I 

I 
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Observando os dados da Tabela 3.14, nota-se que houve 
consistencia. entre os resultados da "proporcao da cultura corrigida" 
com os da "proporcao da cultura obtida atraves de verdade terrestre". 

A baixa percentagcjm de acerto em todos os segmentos 
testes indica que o desempenho do classif icador deve ser melhorado, o 
que provavelmente podera ser feito atraves de outra sistematica, na 
fase de treinamento na qual sao obtidos os parametros para o 
classificador. 

Para fazer uma estimativa da area da cultura da 
cana-de-acucar nos municTpios de Aracatuba e Guararapes teria sido 
necessario que a escolha dos segmentos testes tivesse sido aleatoria 
para poder-se entao agregar as proporcoes obtidas nos segmentos 
testes e, assim, ter-se uma estimativa da area plantada com 
cana-de-acucar nos municipios em estudo. 

Tomando como base a media das proporcoes da cultura 
obtida nos segmentos testes pelo Prooedimento INPE, ter-se-ia uma 
alta precisao na estimativa da area canavieira, na regiao de estudo, 
uma vez que a diferenca relativa entre a media das proporcoes da 
cultura pelo Prooedimento INPK e a media das proporcoes da cultura 
obtida atraves da verdade terrestre (fotografias aereas) foi de 
apenas 2%. 

Embora os resultados de estimativa de area, utilizando 
0 mitodo automatico (Pvoaedimento INPE), tenha produzido excelentes 
resultados neste experimento, deve ser considerado que e a primeira 
vez que tal procedimento foi testado para a cultura da 
cana-de-acucar, nesta regiao, e, portanto, ele seria apenas uma 
primeira aproximacao de uma metcJologia, satifatoria para fins de 
previsao de safras. 




3.4 - MODELOS DE PRODUTIVIDAOE 


3.4.1 - PRODUTIVIDADE ATRAV^S DE MODELO AGROMETEOROLOGICO 

De toda a aroa canavieira das cinco destilarias 
existentes na area piloto, 46,5% era plantada com cana planta de ano 
e meio e 35,5% com cana soca. A area restante (18%) era ocupada com a 
1? e a 2? ressoca e pertencia a duas destilarias apenas. 

Foram obtidos pelo modelo agrometeorologico o 
rendimento maximo (RM) e o rendimento real (RR) para a cana planta de 
ano e meio e para a cana soca, mostrados na Tabela 3.15. Para a cana 
de ano e meio obteve-se a media de RM e RR para o plantio de novembro 
de 1981 a abril de 1982 com urn perTodo de crescimento variando de 16 
a 20 meses. Para a cana soca considerou-se a media de RM e RR para 
cortes da cana entre maio de 1982 e outubro de 1982 com urn perTodo de 
crescimento de 10 a 14 meses. 


TABELA 3.15 


VALORES MEDIOS DE RENDIMENTO MAXIMO (RM) E DE RENDIMENTO REAL (RR) 
REFERENTES A CANA DE ANO E MEIO E A CANA SOCA 



CANA DE ANO-E-MEIO 

CANA SOCA 

Rm VMENTO MAXIMO 

135 ton/ha 

90 ton/ha 

RENDIMENTO REAL 

125 ton/ha 

86,5 ton/ha 
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Na Tabela 3.16, encontram-se os rendimentos medios de 
cana de ano e meio e de cana soca, obtidos junto a cada uma das 

f 

destilarias da irea piloto. 


TABELA 3.16 

RENDIMENTOS MEDIOS EM TON/HA DE CANA DE ANO E MEIO E DE CANA SOCA EM 
CADA UMA DAS CINCO DESTILARIAS DA AREA PILOTO, OBTIDOS ATRAVES 
DAS INFORNACOES DAS DESTILARIAS 



ALCOAZUL 

ARALCO 

CRUZALCOOL 

DESTIVALE 

UNIALCO 

MEDIA 

GERAL 

MEDIA GERAL 
EXCETUANDO 
CRUZALCOOL 

RENDIMENTO MEOIO DE 
CANA DE ANO E MEIO 
TON/HA 

100,0 

130,0 

87,0 

133,0 

112,0 

112,4 

118,7 

RENDIMENTO MEDIO DE 
CANA SOCA 
TON/HA 

100,0 

90,0 

54,0 

75,0 

90,0 

81,8 

88,7 


Os rendimentos de cana de ano e meio e de cana soca, 
obtidos atraves de modelo agrometeorologico e das informagoes 
coletadas nas destilarias, foram calculados utilizando a media 
ponderada dos rendimentos da cana de ano e meio e de cana soca, cujas 
ocorrencias eram de 46,5% e 34,5%, respectivamente. 

A Tabela 3.17, mostra o rendimento geral para a regiao 
de estudo, sem levar em consideragao a cana de 1$ e 2? ressoca. 
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TABELA 3.17 

f 

RENDIMENTOS DE CANA DE AND E MEIO E DE CANA SOCA OBTIDOS ATRAVES 
DO MODELO AGROMETEQROLOGICO E DAS INFORMACOES COLETADAS 
JUNTO As DESTILARIAS DA AREA PILOTO 



VALUR EM TON/ha 

RENDIMENTO MAXIMO FSTIMADO PELO MODELO 

115,5 

RENDIMENTO REAL ESTIMADO PELO MODELO 

108,3 

RENDIMENTO OBTIDO POR INFORMACOES DAS 
CINCO DESTILARIAS 

99,1 

RENDIMENTO OBTIDO POR INFORMACOES DAS 
DESTILARIAS, EXCETO CRUZALCOOL 

105,7 


0 rendimento fornecido pelo modelo superestimou a 
producao em 9,2%, quando comparado com o rendimento obtido por 
infomacao das cinco destilarias. Excetuando o valor de rendimento 
obtido junto a destilaria Cruzalcool, a diferenca entie o rendimento 
obtido pelo modelo e aquele obtido jr to as quatro destilarias cai 
para 2,5%. 


Para o calculo da estimativa do rendimento real 
considerou-se que somente a deficiencia hidrica afetaria o rendimento 
da cultura, nao levando em consideracao, portanto, a infiuencia de 
outros fatores, tais como: tratos culturais, adubacao, ocorrencia de 
pragas e doencas, invarao de ervas daninhas e outros. 

Nos trabalhos de campo observou-se que as tecnicas 
operacionais da cultura de cana nao estavam sendo seguidas com rigor 
e isto, naturalmente, afetou diretamente o rendimento da cultura. 
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Este fato foi mais evidente na destilaria Cruzalcool, razao pela qua! 
0 rendimento obtido por informagao das destilarias e apresentado de 
duas tnaneira's: 

- informagao do rendimento medio das cinco destilarias e 

- informagao do rendimanto medio das destilarias, excetuando a 
destilaria Cruzalcool. 

A quantificagao do efeito dos tratos culturais sobre o 
rendimento da cultura da cai -de-agucar e dificil de ser determinada. 
No entanto, a comparagao entre o rendimenco real, estimado pelo 
modelo, e aquele obtido nas destilarias fornece urn indicative da 
eficiencia da produgao agrTcola da cana-de-agucar na regiao de 
Ara^atuba e Guararapes. 

3.4.2 - '"RODUTIVIDADE ATRAVES DE INDICES VEGFTATIVOS 

Na Tabela 3.18 sao apresentados os coeficientes que 
mostram a correlagao entre os valores de rendimento de cana, dados 
espectrais do LANDSAT e indices vegetativos. 

TABELA 3,18 

COEFICIENTES DE CORRELAgAC (r) ENTRE 0$ VA LORES DE RENDIMENTO D A 
CANA-DE-AgOCAR, DADOS ES PF "RAIS E INDICES VEGETATiVOS 


PASSAGENS 


CANAIS 

DO MSS 


Indices vegetativos 

(IV) 

B 

B 

C 6 

B 

B 

R 

7-5 

IVT 

6 

IVT 

7 

16/Fevereir. , o’ 

r 

- 0,30 

-0,36 

0,28 

-0,42* 

- 0,46* 

0,48* 

0,41* 

0,4^1* 

07/Maio/83 

-0,44* 

-0,46* 


0,002 

- 0,32 

0,51* 

0,43* 

0,51* 

10/Jul ho/83 

-0,29 

-0,25 


0,39 

-0,35 

0,33 

0,30 

0,34 

26/Ju1ho/E3 

-0,40 



0,18 

-0,24 


0,22 

0,27 

27/Agosto/83 

j 

-0,12 



-0,33 


0,43* 

- 0,4^* 


0,43* 
1 

0,43* 






















T 
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Ndo obstante a grande ocorrencia de variedades 

diferentes em um mesmo talhao, pode-se observar, na tabela 3.17, que 
37,5% dos coeficientes obtidos apresentaram-se significativos 
estatisticamente ao nTvel de a = 0,05. Ainda observando a Tabela 

3.17, os indices vegetativos apresentaram, de inodo geral, como 
tnelhores indicadores para estimar o rendiniento de cana-de-acucar do 
que os dados espectrais dos canals do MSS. Quase todos os indices 

vegetativos das passagens de 16 de fevereiro, 7 de maio e 27 de 

^ agosto apresentaram-se significativamcnte correlacionados ao nTvel de 
5%. Isto permite, entao, a utilizagao de indices vegetativos obtidos 
em fevereiro para estimar o rendimento de cana, ja que esta data 
seria muito oportuna para estimar este rendimento pelo menos dois 
meses antes do inicio da saf ra. 

Muito embora os testes de correlagao entre os 
rendimentos de cana e os indices vegetativos tenham-se mostrados, em 
sua maioria, significativos estatisticamente, neste estudo, os 
coeficientes nao foram suf icientemente altos para permitir estimar, 
com precisao, o rendimento da cana em fungao dos indices vegetativos 
obtidos. 

3.4.3 - PRODUTIVIDADE BASEADA NA AREA DE CANA CORTADA 


Os mapas das Figuras 3.4, 3.5 e 3.6 mostram a evolugao 
de areas de cana cortada nos municTpios de Aragatuba e Guararapes 
reterentes aos meses de julho, agosto e setembro, obtidas da 
interpreta;ao de imagens do LANDSAT. 

Areas de cana cortada apresentam-se em tonal idades 
Claras em imagens do canal 5 do sensor MSS/LANDSAT, em contraste cc 
a cana nao-cortada que se mostra en^ tonal idades de cinza. 





Fig. 3.5 - Areas de cana cortada ate 27 de agosto de 1983 
nos municTpios de Aragatuba e Guararapes, obtJ_ 

das da interpretagao de imagens do LANDSAT. j 


i 



ORIGINAL PACI sS 
OF POOR QUALITY 



Fig. 3.6 - Areas de cana cortada ate 2S de setembro de 1983 
nos municTpios de Ara?atuba e Guararapes, obti 
das da interpretagao de imagens do LANDSAT. 
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A label a 3.19 apresenta os resultados de area de cana 
cortada obtidos de- satelite e de aeronave, tonelagem de cana moida 
pelas destilarias e o rendimento medio baseado na area de cana cortada 
das cinco destilarias em funcionamento nos dois municTpios. 

TABELA 3.19 

RENDIMENTO MEDIO DE CANA DAS DESTILARIAS DE ARACATUBA E GUARARAPES 


DATAS 

Area de cana 
CORTADA (ha) 

producao 

RENDIMENTO MEDIO 
(T/ha) 

22.06.83 

2.167* 

200.074 

92,3 

10.07.83 

3.668 

350.991 

95,7 

26.07.83 

4.723 

469.238 

99,3 

27.08.83 

8.157 

808.949 

99,2 

28.09.83 

13.319 

1.112.062 

83,5 


* obtida atraves de fotografias aereas. 

Os menores rendimentos alcancados ate junho podem ser 
atribuTdos ao uso de cana cortada para o p'lantio de outras areas de 
cana. Pela analise do calendario agrTcola desses municTpios 
observou-se que 80% do plantio de cana em 1983 ocorre nos meses de 
abril , maio e junho. 

Poster iormente, os rendimentos medios calculados no 
final de julho e agosto mantiveram-se quase constantes, voltando a 
declinar quando se trataram informacoes do final de setembro, 
apresentando uma diferenqa relativa de 15,8% em relacao ao mes 
anterior. Possivelmente fatores agronomicos tenham influenciado na 
queda de rendimento da cana-de-acucar que foi cortada durante o mis de 
setembro. 



CAPTTULO 4 
CONCLUSOES 


4.1 - AO NIVEL MUNICIPAL 

a) Tra^ado da rede hidrografica e rodoviaria 

- Com relagao ao tracado da rede hidrografica verificou-se 
nao haver dificuldade na sua identif icacao atraves de 
imagens de satelite. 

- Quanto ao tracado da rede rodoviaria as imagens MSS do 
LANDSAT apresentaram limitaqoes na sua identificacao, 
principalmente para as vias de acesso secundario. 

b) Mapeamento e estimative da area de classes de uso do solo 

- As imagens LANDSAT apresentaram-se viaveis para o 
mapeamento e a estimative da area de classes de uso atual 
do solo. 

c) Desenvolvimento de modelos para estimative de produtividade 

visando previsao de safras 

- 0 modelo agrometeorologico forneceu urn indicative do 
desemptnho do si sterna de producao da cana-de-acucar na 
regiao de estudo e resultou em dados de produtividade 
bastante proximos dos estimados pelas destilarias. 

- Os dados MSS do LANDSAT, transformados em indices 
vegetativos, apresentaram boa correlapao com os dados de 
campo de Liomassa da cultura, e que indica potencial de 
utilizacao em futures modelos de produtividade. 

- 0 uso de imagens do LANDSAT para avuliacao de area de cana 
cortada, aliado a informaca> de quantidade de cana 
recebida pela destilaria, permitiu a estimativa do 
rendimento da cultura de forma satisfatoria. 
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4.2 - AO NTVEL DE PROPRIEDADE 

a) Viabilidade de localizagao e delimita?ao de propriedades 
agrTcolas nas imagens do LANDSAT 

- Dentre os diversos produtos fotograficos e topograficos 
utilizados na delimitagao de propriedades agrTcolas, 
concluiu-se que o melhor produto para cadastramento foram as 
fotografias aereas de arquivo (1978 e 1972). 

- Atraves de transposigao dos limites das propriedades 

agrTcolas para as imagens do LANDSAT foi possTvel local izar 
e identificar as propriedades de forma satisfatoria em 
imagens de satelite. 

b) Estimativa da area da cana-de-agucar 

- A percentagem de cana interpretada corretamente atraves de 
imagens do LANDSAT foi de 63,12%. A utilizagao de uma 
abordagem alternativa que envoi ve tanto o uso de imagens do 
LANDSAT quanto o trabalho de campo elevou este Tndice para 
91,77%. 

c) Com relagao ao tamanho mTnimo de propriedade, nao houve 
dificuldade com relagao a sua identificagao em imagens 
do LANDSAT na escala de 1:250.000 sendo que a menor propried^ 
de com financiamento para implantagao da cultura canavieira 
existente na regiao de estudo foi de 40ha, enquanto a area mT 
nima de cana detectavel foi da ordem de 25ha. Entretanto, d^ 
pendendo do contexto, talhoes menores podem ser perfeitamente 
caracterizados. 


i 




- 85 - 


CAPTTULO 5 
CONSIDERACOES FINAIS 


Muita experiencia foi adquirida durante o 
desenvolvimento deste projeto, principalmente devido a abordagem de 
varies problemas inerentes a fiscal izacao das operacoes do credito 
agrTcola do Banco do Brasil S.A. 

0 cadastramento de propriedades agricolas, praticamente 
inexistente para o Estado de Sao Paulo, foi urn dos primeiros pontos a 
ser estudado. Chegoi-se a conclusao que o simples uso de fotografias 
aereas convencionais de arquivo e trabalhos expedites de campo podem 
fornecer os limites dos imoveis agricolas com boa precisao. A 
transposicao destes limites para imagens na escala de 1:250.000 tambem 
se mostrou adequada. Desta maneira, seria possivel ter, com urn custo 
relativamente baixo, o cadastre em forma cartografica de propriedades 
agricolas. Urn treinamento rapido de fiscais ou auxiliares das agencias 
do Banco do Brasil S.A. sobre o manuseio de fotografias aereas 
convencionais seria suficiente para torna-los aptos a levantar os 
limites de imoveis rurais junto aos seus proprietaries por ocasiao do 
pedido do financiamento. 

No tocante a estimativa de area da cultura canavieira 
utilizando imagens do LANDSAT, encontraram-se muitas dificuldades para 
sua identificaqao. Isto sugeriu a adocao de novas estrategias 
metodologicas , nunca testadas anteriormente , como o uso exaustivo de 
dados historicos do LANDSAT, desde o inicio de implantapao da cultura 
na area piloto (1980) ate 1983, auxiliado por varias informaedes de 
campo. Isto permitiu que se alcangassem resultados bastante 
signif icativos de avaliagao de area do tema cana. 

Com relacao a produtividade, foram testadas tres 
abordagens para a cultura da cana-de-agucar. Os resultados, 
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considerados positives, enfatizam a continuagao das pesquisas no 
intuito de controlar melhor as variaveis estudadas, incluir outras e 
testar os procedimentos em diferentes regioes. 

A tecnologia de sensoriamento remote tern avangado 
rapidamente nos ultimos anos e, seguramente, a comunidade de usuarios 
de dados de satelites de recursos terrestres tera informagoes 
atualizadas com alta frequencia de imageamento por muitos anos. 

Neste projeto, util izaram>se dados de sensor MSS/LANDSAT 
com uma resolugao no terrene da ordem de 0,48ha. 

0 quinto satelite da serie LANDSAT ja esta em orbita 
desde margo de 1984, levando a bordo outro sensor - o Thematic Mapper 
(TM), com uma resolugao de 0,09ha e operando em sete bandas espectrais 
(4 no MSS). 


0 satelite frances denominado SPOT (Sisteme Probatoire 
d 'Observation de la Terre), com previsao de langamento no primeiro 
semestre de 1985, tera dois modes de operagao: urn, multiespectral de 
tres canais, de resolugao de 20m e o outro, pancromatico (apenas urn 
canal) de resolugao de 10m. 

Pode-se afirmar a prio*"! que este aumento de resolugao, 
tanto espacial quando radioiuitrica, vai beneficiar, em muito, a 
agricultura no que diz respeito as previsoes de safras agricolas, 
tanto na identificagao e avaliagao de areas de varias culturas quanto 
na obtengao de dados espectrais mais precis'^s que alimentarao os 
model os de produtividade em desenvolvimento. 

Este projeto, embora de pequena duragao, gerou 
resultados que efetivamente encorajam a continuagao de esforgo de 
pesquisa com o objetivo de aprimorar os metodos descnvolvidos. 
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A disponibil idade de sistemas sensores de maior poder 
resolutivo .(LANDSAT-5 TM), aliado a um processo de cadastramento de 
propriedades com a indicacao de classes de uso da terra por ocasiao do 
processo de financiamento, e a experiencia adquirida neste projeto 
permitem antecipar que os resultados poderao ser sensivelmente 
melhorados. 


A continuacao desse trabalho devera envoi ver a regiao 
piloto estudada (Aracatuba e Guararapes) por ja possuir como base de 
apoio substancial para desenvolvimentos metodolSgicos. Paralelamente, 
um esforco de transicao de uma fase experimental para uma operacional 
deveria ser colocado envolvendo municTpios vizinhos. 
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